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ABSTRACT 


A Combined Helmholtz Integral Equation Formulation (CHIEF)-defined 
enclosing sphere placed around an acoustic projector is used to calculate far field 
response data from near field measurements. Pressure response data at a specified 
frequency is obtained from a near field linear array. Helmholtz integral relations for the 
enclosing sphere and integrals of the free-space Green’s function and its gradient for 
defined near field point locations are used along with physical assumptions to form an 
overdetermined system. The overdetermined system is solved via least squares yielding 
values of pressure and velocity corresponding to defined locations on the enclosing 
sphere. The enclosing sphere’s values of pressures and velocities are then used with 
integrals of the free-space Green’s function and its gradient to calculate far field 
response. 
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I. BACKGROUND 


The testing and eharaeterization of U.S. Navy underwater acoustie devices 
remains a fundamentally important activity linked directly to both fleet and mission 
readiness. Testing of new prototypes, sampling of production line devices and re¬ 
qualification of in-service assets are all components of this readiness. These testing 
efforts represent measures put in place to help assure that the required level of acoustic 
performance is being achieved. 

The Underwater Sound Reference Division (USRD) represents the Navy's 
primary activity for underwater acoustic calibration, test, and evaluation measurements. 
In this capacity, the USRD maintains specialized measurement facilities with the 
capability to simulate real-world ocean environments. The USRD is part of the 
Submarine Sonar Department of the Naval Undersea Warfare Center Division Newport 
located in Newport, Rhode Island ( www.npt.nuwc.navy.mil/USRD) .^*^ 

Within the USRD is an acoustic measurement facility designated as the Acoustic 
Pressure Ta nk Facility (APTF). The APTF test vessel consists of a closed steel tank 3.81 
meters in diameter and 11.1 meters in length with two access ports. The facility allows 
the acoustic test and evaluation of devices and materials over a wide range of simulated 
ocean environments. Temperatures from 2 degrees Celsius to 35 degrees Celsius and 
hydrostatic pressures up to 18.6 megapascal (a depth equivalent of 1860 meters) can be 
achieved. 



Naval Undersea Warfare Center Division Newport - Aconstic Pressnre Tank Facility 
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The APTF represents a significant U.S. Navy capital investment. As Navy 
demands and needs have changed over the years, it has become increasingly more 
important to revitalize existing Navy assets to meet new requirements. While it is a 
common need to characterize the performance of acoustic devices under realistic ocean 
environmental conditions, in some instances these devices may require test distances in 
excess of that achievable within the confines of existing Navy test vessels. As an 
example, these devices may entail designs that have large apertures yet may also function 
at relatively high frequencies (as the U.S. Navy High Frequency Sail Array). This 
distance required for proper acoustic characterization is designated as that distance where 
far field conditions occur. Test distances that meet or exceed this far field criterion 
allow the device to function as it would in the free-field condition of the open ocean. 
Testing at distances that fall short of this criterion are designated as being tested in the 
near field. Testing in the near field is subject to contaminations that obscure the true far 
field performance of the device and thus is typically not desirable. 

This thesis effort addresses a technique that uses near field acoustic data with a 
suitable math model to extrapolate far field performance. The use of “CHIEF for 
MATLAB”'-^^ coupled with a least squares solution methodology are utilized in this 
technique. The data output rendered from the technique mimics data expected within test 
vessels that would be prohibitively expensive to fabricate and maintain. 

The remainder of this thesis is organized as follows: 

Chapter II addresses the definition and meaning of near-field and far- field conditions. 

Chapter III provides an overview of the Combined Helmholtz Integral Equation 
Eormulation (CHIEE) and the A and B matrices that are used in subsequent computations 
in the thesis methodology. 

Chapter IV details the steps of the thesis near-field to far-field methodology and criteria 
for assessing the data output. 
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Chapter V provides an example solution set using the methodology. 


Chapter VI provides an error analysis of the technique when positional uncertainties of 
the various acoustic devices used during a test exist. 

Chapter VII looks at other specific test examples inclusive of testing in more extreme 
near-field cases and test results expected from what is likely to be considered as a typical 
test set-up and scenario. 
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II, NEAR-FIELD/FAR-FIELD OVERVIEW 


A, INTRODUCTION 

A fundamental physical dimension addressed in conducting acoustic 
measurements is separation distance. The desired separation distance between devices in 
a test setup places them at distances where 1/r spreading occurs (with r defined as that 
distance between the projecting and receiving devices). This mimics expected 
performance in a free-field state and enables an accurate characterization of the device 
for the typical ocean environment for which it is designed. 1/r spreading gives the 
expected 6 dB loss when separation distance is doubled This type of spreading occurs 
when “far-field” conditions have been met. 

Testing of a transducer with finite dimensions (i.e., not a point source) at 
distances inside the far-field limit results in performance characterization that is obscured 
by phase superposition of the comprising elements. That is, the complex response 
measured at a particular point is made up of a number of individual responses each with 
widely varying phase terms. Thus, the “near field” of a particular set-up is characterized 
by radical fluctuations in pressure with changes in range (in a given direction). This 
obscures the nature of the device’s true far field performance. 

The following curve demonstrates the variation in the near field and far field 
responses seen for a particular acoustic device (a linear array). The response curve was 
obtained by summing up the complex responses from all source elements at each of 
several specified ranges. Note the widely changing response in the near field of the 
device and the 1/r spreading resulting at the farther ranges. 
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Figure 1 - Pressure Response with Change in Range 


The field data at varying ranges in the previous figure (Figure 1) was computed 
for a 30.5 cm linear array using a test frequency of 600 kHz. A far field distance of 
approximately 37 meters is estimated for this device. The data set was normalized to the 
maximum response computed. Computations were based upon modeling the source as a 
“Continuous Line Source” as detailed in Kinsler and Frey, Section 7.3 (4'*' Edition) 

B. THEORY 

The properties of the corresponding acoustic field for a given device can typically 
be estimated from the geometry of the source and its wavelength, X , of operation. The 
extent of the near field can be estimated through the examination of the acoustic 
contributions from the device’s various geometric components. The figure that follows is 
used to examine these details. 
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Figure 2 - Geometry used in estimating extent of Near Field 


In the above figure, , represents the distanee from the furthest element of the souree 

(in this example this is a line array of length L) to a field point. The magnitude, , 
represents the distanee from the nearest element of the souree to the same field point. 
The differenee between these two magnitudes is defined as A r = . At great 

distanees, this differenee (Ar ) is minimal. At eloser distanees and as the field point 

approaehes the souree on axis, this differenee ean be signifieant. When the differenee 
approaehes XU (over the differenee obtained at great range), the phases of the signals 
from the individual points on the souree eombining at the field point will be of 
suffieiently different value. The net effeet of this is that the axial pressure. Pax, can 
ehange radieally over small ehanges in distanee (see Figure 1).'-"^’^*’^ 

Referencing Figure 2 again gives the following expression for indicating a 
reasonable boundary for expected far field performance: 

r_/L = L/4A (1) 
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L/2 represents half of the maximum dimension of the souree. 


A rule of thumb for far field separation distanee that has been used in the USRD faeilities 
has been based upon 

D " 

( 2 ) 

with D equal to the maximum dimension of the aperture 
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Ill, CHIEF - INTRODUCTION 


A, BACKGROUND 

The CHIEF (Combined Helmholtz Integral Equation Formulation) program was 
originally developed at the Naval Ocean Systems Center (NOSC) from 1965-1970 by 
Drs. H.A. Schenck and G.W. Benthien. It was originally developed to compute acoustic 
radiation that would be produced from a given arbitrarily defined body. Over the years 
its capabilities have been expanded by others to utilize more powerful and available 
computing hardware and to include plane-wave scattering and finite-element defined 
input descriptions . 

The version of CHIEF used in this effort was written for implementation with 
MATEAB by Mr. Stephen Forsythe, NUWCDIVNPT, during fiscal years 1999-2000. 
This represents a unique implementation of CHIEF that takes advantage of the matrix 
handling capabilities of a commercial off the shelf product. While achieving the 
functionality of the original CHIEF program, Mr. Forsythe further augmented this 
version to utilize spatially defined material properties rather than just relying on boundary 
value problem descriptors . The open nature of the coding that was implemented also 
allows access to non-traditional CHIEF outputs. It is this latter implementation detail that 
made this version of CHIEF the proper choice for this effort. Specifically, the extraction 
of the CHIEF A and B matrices for use outside of the protocols of CHIEF was the unique 
feature that this version provided and that this thesis methodology required. 

B, CHIEF THEORY 

This section provides a brief glimpse at the definition and formulation of the 
CHIEF A and B matrices and the use of the Helmholtz Integral formula to deduce 
pressure at a point in space. This thesis effort does not concentrate on the exacting 
derivation of these parameters but rather provides an overview. More complete details 
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may be found on this derivation in the doeumentation provided with the CHIEF 1988 
release . 


For both aeoustie radiation and aeoustie seattering problems, Green’s theorem 
speeifies that the total aeoustie pressure p ( pinc + Ps, with pinc defined as the ineident 
pressure on a surfaee S due to sourees in the region of interest and ps defined as the 
eomplex seattered pressure exterior to S) and the normal veloeity v satisfy the Helmholtz 
integral formulas: 



p{x)-p inc (x) ; X exterior to 5* 
^ p{x) - pinc{x) ;xonS 
— pinc (jt); X interior to S 


with: 


1 

G(x,x’) = --^ 

4;r x-x 



( 2 ) 


where: 

G(x,x') = three dimensional free-spaee Green’s funetion and 

k: aeoustie wave number 
co: radian frequeney 
e: sound speed in the fluid 

10 



-: signifies differentiation in the outward normal direetion 

dn 


and a harmonie time dependenee of e'® is assumed. 

If ^ is some point loeated on the surfaee, equation (1) from above leads to: 

Picr)^G{C,^)dsi'^) = Pinc{C) + i(^p\ v{cr)G{(;,cr)ds{<j) (3) 

C. CHIEF - NUMERICAL APPROACH 

Sinee the goal in CHIEF is to solve this problem numerieally, the surfaee S is 
broken up into an Ns number of non-overlapping “patehes”. The pressure and veloeity at 
eaeh of these patehes is approximated by average values, designated as Pn and Vn, 
respeetively. 

The diseretization of equation (3) leads to: 

1 v, d 

tr- S^. or)ds{cj) = ) 

2 „=i on^ 

+ fJ G(C , cj)ds{cj) (4) 

«=i s„ 


The matriees A and B are subsequently defined as follows: 
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A 


mn 


2 





d 

dn^ 




B 


mn 



G{<^„,a)ds{a), 


(5) 


If the variables p , v, and pi„c are defined as vectors of surface pressures, 
velocities and incident surface pressures, then equation (4) can be more simply stated as: 

A p = B V + Pine (6) 

Note: further refinement of the A and B matrix definitions in equation (5) above 
are required to handle specific wave numbers and surface pressure conditions. The 
matrix A is otherwise not suited to achieving an over determined system solvable via 
least squares solution technique. This further refinement is not addressed here. 

As mentioned in the early part of this section, the Helmholtz Integral relations can 
also be used to deduce pressures at defined locations in space. The pressure at any field 
point X (near field or far field) can be expressed in terms of the surface pressure and 
normal velocity per the following relationship: 




(V) 



G(x, a) + i(i>pv{<j)G{x, a) \ds{<j) 


In like fashion as the formulation of the CHIEF A and B matrices, it desired to 
solve equation (7) numerically by breaking the surface S into a mosaic of non¬ 
overlapping patches. This leads to the expression identified in equation 8 that follows. 
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Psi^) = Pi^)- 


( 8 ) 


" r d " r 

= — G{x,CT)ds{CT) + iQ)p'Yv„\G{x,CT)ds{cT) 

n=l n=l g^ 

Using the condition where pmc{x) is equal to zero (indicating that there is no 
incident pressure on the surface), the discrete notation of equation 8 may be further 
simplified to: 

p{x) = p delG + V G (9) 

where each of the constituents (from left to right) on the rhs of equation 9 are vectors of 
the pressure at a known surface location, the integration of the normal derivative of the 
Greens function for the given set-up, the normal velocity at a known surface location and 
the integration of the Greens function for the given set-up, respectively. Please note that 
X designates the three dimensional location in space of a defined field point and cr 
(equation 8) designates the three-dimensional location in space of the surface patches. 

D, SUMMARY 

CHIEF utilizes the Helmholtz integral relations of surface pressures and normal 
velocities along with a defined surface incident pressure condition. Computations within 
CHIEF proceed numerically via the re-structuring of the continuous surface into a mosaic 
of individual non-overlapping surface patches. This re-structuring allows the formulation 
of the generalized equation: 

Ap = Bv + pinc (10) 
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with p , V, and pinc defined as veetors of surfaee pressures, normal veloeities and ineident 
surfaee pressures, respeetively. A is the matrix formed from surface integrals of the 
normal derivative of pressure over the surface patches and B is the matrix formed from 
surface integrals of the pressure over the surface patches. Both the A and B matrices can 
be formed from known acoustic constants and a defined geometry . 

CHIEF also yields a mechanism to calculate pressure at a defined field point via 
the use of known surface pressures and normal velocities and integrated three- 
dimensional Greens function derived relationships (between surface locations and the 
defined field point x). This leads to the following relationship: 

;7(x) = p delG + V G (11) 

with p{x) defined as the pressure at a designated field point x, p and v defined as vectors 
of surface pressures and normal velocities, respectively, and delG and G defined as the 
integrated Greens functions relationships between field point x and the surface patch 
locations. 

Equations 10 and 11 in the above summary form most of the basis that will be 
utilized in setting up the near field-far field extrapolation problem’s defining equations. 
The methodology will require that Equation 10 be equated to zero (indicating that no 
external sources are present) and equation (11) will be used to relate to pressures at 
defined near field points. 

Due to the numerical approach that CHIEF utilizes, it is important that the number 
of patches used in the mosaic to cover the surface, S, be properly sized. A single 
complex value of pressure and velocity per patch location is a required condition in 
CHIEF; these values of pressure and velocity over the entire patch are represented by 
values computed for the patch center. 
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IV. STEPS OF THE METHODOLOGY 


A. INTRODUCTION 

This section outlines and details the steps of the methodology required to take 
near field array element data through the processing steps required to achieve an 
extrapolated far field response. The steps provide the procedural details as well as criteria 
for parameter selection. Practical considerations of limitations of the processes are also a 
necessary part of these steps since the effort entails the implementation of the technique 
on available computing platforms. This section is supplemented by the “Explanation of 
MATLAB Coding” section that appears in Appendix C. 


CHIEF Enclosing Sphere placed around Linear Transmitting Array 


BLACK: Transmitting Array 

RED: CHIEF Enclosing Sphere .. 
BLUE: Reicjeive Array 



y (meters) 


-2 -2 


X (meters) 


Figure 3 - Typical Methodology Set-up showing device under test, receive array and 
CHIEF Enclosing Sphere 
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The figure above shows a typical set-up that forms the basis of the methodology. 
The transmitting array and the receive array both represent actual physical devices, while 
the CHIEF enclosing sphere represents an artificial surface in space. The methodology 
seeks to define the pressure and velocity characteristics of this spherical surface for 
eventual extrapolation to far field pressure. 

While some of the early steps are not execution order specific, it is practical to 
treat all of these steps as sequential in the order that they are presented. 

B, DEFINE CHIEF ENCLOSING SPHERE 

This step of the process takes into account the physical size of the acoustic device, 
and its frequencies of operation. Since the enclosing sphere (formulated using “CHIEF 
for MATEAB” protocols ) will be formed from multiple repeated copies of a single 
“slice” of a sphere (imagine a segment from a section of a tangerine), an estimate of the 
number of comprising slices required, ie: its symmetry order (sph sym), is performed in 
this step. The number of divisions within thi s single “slice” of the sphere (designated as 
patches) are also formed in this step. The slice will, therefore, be represented by a 
sph_numv x 1 matrix with sph_numv representing the number of equally spaced patches 
in the direction perpendicular to the symmetry order . (Eimitations are imposed on both 
of these quantities based upon computational capabilities and efficiencies of both the 
MATEAB coding and the standard COTS computer hardware used for this effort.) 

C. DEFINE RECEIVE ARRAY AND PROXIMITY TO ACOUSTIC DEVICE 

This step defines the number of receive array elements (numb) and the array’s 
distance from the acoustic device under test. It is assumed that the elements in this array 
are in a vertical line and are equally spaced. The technique does not necessarily require 
this for any particular technical reason, rather this is based more upon the expected 
configuration that will be suited for an experimental set-up. Maximum spacing of the 
elements in this array is expected to be dependent upon the required frequencies of 
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operation. The maximum number of receive elements for this technique will be 
established more by logistical considerations rather than for technical reasons (ie: there 
will be some practical limit to the number of elements that can be made available for the 
receive array based upon available sensors in inventory and the number of channels 
available on the data acquisition system). The distance of the receive array from the 
acoustic device will be in the near field of that device. Part of the evaluation of this 
technique will be to examine final extrapolated data to assist in determining the 
robustness of the technique with decreasing separation distance from the acoustic source 
and varying receive array element spacings and numbers. This technique will require that 
the receive array be placed outside the boundaries imposed by the enclosing CHIEF 
sphere. 

D, RECEIVE ARRAY MEASUREMENTS/SIMULATION 

This step involves the acquisition (or simulation) of complex receive response at 
each of the receive array elements. This step requires that data be collected for each of 
the receive array elements at acoustic device rotation angle intervals corresponding to 
(360 degrees /symmetry order of the enclosing sphere). The sampling (or later 
interpolation) of this data at the prescribed angles allows a direct relationship between the 
data collected and the defined “slices” of the enclosing sphere (established by the 
selected symmetry order). 

E. COMPUTE ORIGINAL A AND B MATRICES VIA CHIEF FOR MATLAB 

This step involves the computation of the original A and B matrices using 
“CHIEF for MATLAB” tools . The A and B matrices are formed from the surface 
integrals in CHIEF. The following is taken from documentation provided for “CHIEF for 
MATLAB” by S.E. Forsythe (NUWCDIVNPT): “The standard problem statement for 
any CHIEF problem is: Ap = Bv + pinc, where A is the matrix formed from surface 
integrals of the normal derivative of pressure over the surface patches, B is the matrix 
formed from surface integrals of the pressure over the surface patches, pinc is the vector 
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of incident pressures on the surface patches due to sources in the region of interest, p is 
the (to be solved for) vector of total pressures on the surface, and v is the (to be solved 
for) veetor of veloeities on the surface.” 

CHIEF for MATLAB has the capability to provide the A and B matrices and 
eease further eomputations (to solve for resulting pressures and veloeities). This is 
desired for our far field extrapolation teehnique since the A and B matrices will be 
augmented with other faetors in the next step. The primary inputs to ereate the A and B 
matrices depend upon the aeoustie parameters of sound speed, density of the test medium 
and frequency of interest and upon the geometric configuration of the surface of interest. 
In our case, the surface of interest is the enelosing sphere from the first step. Thus, given 
a defined enclosing sphere, the matriees A and B can be formed. 

F. COMPUTE AUGMENTED MATRICES Aa , Ba AND Ma 

This next step begins the departure from the usual manner in whieh the CHIEF 
program is utilized. That is, typieal CHIEF applieations require the definition of a 
transdueer geometry followed by the characterization of a pressure, veloeity or 
impedance across all defined patches. Specification of one of these parameters (pressure, 
veloeity or impedance) is required to make a reduced system (from the original defining 
equations of Ap = Bv + pinc that is uniquely solvable. 

In this thesis methodology, the only eompletely defined geometries are 
represented by the CHIEF enclosing sphere (an artificial surface) and the physieal 
loeation of the reeeive array elements; the geometry of the device under test is not 
required to be known exactly (with the exception that its largest dimension must be 
known). This flexibility to work with a minimal set of known transdueer geometries was 
a desired feature of the methodology. 
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The other signifieant differenee in this CHIEF usage is that the pressure and 
veloeity terms at eaeh pateh loeation on the defined CHIEF enelosing sphere are both 
unknown values that must be solved for using additional information eonstraints. 

A least squares solution is the ultimate goal of the methodology and relies upon 
the use of field values of pressure as the starting point. The relation between the CHIEF 
A and B matriees and pressure and veloeity terms form the fundamental relationship of 
0 = -A p + B V and is examined first below. Since this system by itself is 
underdetermined, the system is augmented by other relationships to achieve 
overdetermined status. The augmentation relies upon the CHIEF A and B matrices and 
their relationship to field pressures; reliance is also made upon expected adjacency 
conditions of pressure and velocity between patches. 

The A and B matrices formed in the last step provide the starting point to develop 
the following relationship: 

Ma = -AaP + BaV 

where: 

Aa = Augmented A Matrix 
Ba = Augmented B Matrix 

p = pressures at all “patches” on the enclosing sphere (a vector) 

V = velocities at all “patches” on the enclosing sphere (a vector) 

Ma = Augmented Field Points Matrix (complex responses at defined points in the test 
medium) 

The augmentation of the matrices in this step puts the measured data and defining 
physical equations in a form that will be suited for a later least squares solution 
determination. The augmentation of the A and B matrices is examined first and is 
described as follows: 
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1) 0 = -A p + B V 


This relationship states that the eomputed A and B matriees (from the geometrieal 
and aeoustie eonsiderations of the previous step) when multiplied by their respeetive 
pressures and veloeities result in zero net pressure on the surfaee when there is no souree 
in the region of interest as required by the Helmholtz theorem. 

2) pressure at a field point = delG p + G v 

This defining equation relates the Greens funetion and its gradient along with 
veloeity and pressure values at eaeh physieal location on the enclosing sphere to complex 
response values measured at (near) field test points. The Green’s function and its 
gradient are determined via geometrical considerations (ie: the physical location of the 
field points with respect to the enclosing sphere). The pressures at field points are as 
measured (or simulated) by the receive array. As in the first step of this augmentation, 
this defining equation utilizes the (to be solved for) pressures, p, and velocities, v, along 
with factors that can be determined ( G and delG). 

3) pressure/velocity equivalence at adjacent patches 

This relationship implies that due to the relatively small size of the patches 
comprising the surface of the enclosing sphere, it is reasonable to assume that adjacent 
patches do not differ significantly in pressure or velocity from one another. They can, in 
fact, be considered to be approximately equal. This allows a physical attribute to be 
associated with the (to be solved for) pressures and velocities. Enforcing this condition 
with a small weight in the least squares solution prevents the extreme case where non- 
realizable pressure and velocity swings are required for mathematical realization 
Thus, this step essentially states that pi = pi+i and Vi = Vi+i, where i represents some patch 
location and i+1 represents the patch next to it in the vertical direction. 
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The resulting augmentation of the A and B matrices, therefore, yields the 
following constituents: 


Aa Matrix 

Ba Matrix 

-Ap 

Bv 

delG p 

Gv 

Pi = Pi+i 

Vi = Vi+I 


Note that pressure terms are related to the Aa matrix while velocity terms are 
related to the Ba matrix. (The details of the actual construction protocols of these 
matrices are examined later.) 

The construction of the augmented field point array, Ma, is now examined. The 
Ma matrix represents the expected outputs or measured values corresponding to the 
constituents of the augmented Aa and Ba matrices. For the three defined parameters 
discussed above, Ma takes on the following values: 

1) Ma= 0 

for the (-Ap+Bv) terms 

2) Ma = complex response measured (simulated) at each element of the receive 

Array for the (delG p + G v) terms 

3) Ma = 0 

For the adjacent equivalence terms of pressure and velocity represented by 
(pi - pi+i) and (vi - Vi+i) 
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Since a least squares solution teehnique will be utilized in the next step, an 
advantage in optimizing the solutionis gained by sealing the matrix constituents. These 
re-scaled matrix components still aeeurately represent the conditions and equalities of the 
set-up. They do, however, de-emphasize the wide variation in magnitudes of the various 
quantities and allow the solution criteria to proeeed based upon more equal weightings. 
This re-sealing is also appropriate to identify as putting the defining equations in 
dimensionless form. 

For example, the Green’s Funetion integral, G, for a partieular set-up returns an 
absolute maximum value of 573 while the gradient returns an absolute maximum value of 
3.82e-4. Adjusting the Green’s Function integral elements by first dividing all G 
elements by rhoc (1500000) and then dividing all matrix elements of (the new) G and 
delG by 3.82e-4 gives matrix elements on both G and delG that have maximum values of 
approximately one. 

While the above protoeol establishes more or less equal weightings, the 
methodology does allow the user the opportunity to place specific emphasis on the 
contributions of the theoretically defined Helmholtz integral relations or of the 
empirically acquired near field measurement data. This is embodied in the alpha and beta 
terms defined below. 

For this effort, the comprising matrices are scaled by the following faetors: 


Matrix 

Scaling Factor 

A 

alpha 

B 

alpha/rhoc 

delG 

beta 

G 

beta/rhoe 

M 

beta 
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with: 

alpha = alphaO / (maximum absolute value in A) 
beta = betaO / (maximum absolute value in delG) 


and: 

alphaO and betaO are as defined by the user to emphasize the theoretieal terms 
(Helmholtz integrals relations) or emphasize the empirieal data (near field 
measurements), respeetively, in the final solution. 

This eompletes the augmentation of the defining matriees. It is signifieant to note 
at this point the expeeted size of the augmented matrix eonstituents. This information is 
tallied below: 

Matrix Size (rows x eolumns) 

A numv x (sph_numv * sph_sym) 

B “ “ 

delG numh x (sph_numv * sph sym) 


Pressure equivalenee ( sph numv-l) x (sph numv * sph sym) 

Veloeity equivalenee “ “ 

Ma ( sph_numv + numh + 2*sph_numv -2) x 1 


The matrix sizes presented here do not eonstitute the full matrix, rather they are 
representative of one “sliee” of the full eirculant matrix . The full matrix size is larger 
by a faetor of sph sym for eaeh of the above eonstituents. 
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The full Aa and Ba matrices are, therefore, of size; 


(sph_sym * ( numh + 3*sph_numv -2) ) x (sph_numv * sph_sym) 

while the full Ma matrix is of size: 

( sph_sym * ( numh + 3*sph_numv -2) ) x 1 

G. SOLVE FOR PRESSURES AND VELOCITIES VIA LEAST SQUARES 

This step uses the least squares solution technique to compute the pressure and 
velocity values at each of the enclosing sphere patches. The output for this step consists 
of (numv * sph sym) pressure values and the same number of velocities. 

The least squares technique allows an overdetermined system (“m” equations in 
“n” unknowns, with m>n) to be evaluated in an optimized manner. This optimized 
manner results in values of x that minimize the quantity |Ax - b|^ for the originating 
equation of Ax = b. 

The vector x that minimizes |Ax - b| is the solution to the normal equations: 

A^ Ax = A^ b 

This vector x = (A^ A)‘^ A^ b is the least squares solution to Ax =b. 

The implementation of the least squares technique for this thesis effort substitutes 
a matrix that will be designated as the AB_Matrix for A and the full Ma matrix for b. 

The AB_matrix consists of the full matrices Aa and Ba that have been interleaved as 
follows: 

Aai Bai Aa 2 Ba 2 Aa^ Ba 3 ... Aan Ban; with n = (sph_numv * sph_sym) 
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The AB_Matrix is actually a circulant matrix of size: 

(sph_sym * ( numh + 3*sph_numv -2) ) x (2 * sph_numv * sph_sym) 

The justification for the interleaving of the Aa and Ba components will be examined in 
detail later. 

The least squares solution, x, given both AB_Matrix and Ma , is a vector of 
alternating pressures and velocities (of total length: 2 * sph_numv * sph sym). These 
values correspond to the pressures and velocities at each patch on the enclosing sphere. 

H. COMPUTE Ma-new MATRIX 

This step takes the pressures and velocities computed in the last step and 
substitutes them into the following equation to compute the new augmented field point 
matrix, Ma-new • 

M^a-new Aa p + Ba V 

The Aa and Ba matrices are the same as that defined previously in step #5. The p 
and the v represent the pressure and velocities computed via the last step. 

This newly computed matrix, Ma-new , can then be used to evaluate the validity of 
the enclosing sphere pressures and velocities, p and v, and whether the solution is likely 
to be properly optimized. This is accomplished by comparing Ma-new quantitatively to 
the original Ma matrix. 

Since Ma-new Utilized a data set derived from a least squares solution (the p and v 
values), it is not expected that it will match up exactly one for one with the originating 
Ma matrix. There are important considerations to be made, however. 
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Before listing this eriteria, a short summary of the eonstituents of the Ma matrix is 
first reiterated. The Ma matrix is divided up into four seetions. 

1) The first section has values corresponding to the (-Ap + Bv) calculations. 
Ideally these values should be zero. 

2) The second section contains the complex response data at each element of the 
receive array (these are near field data points). These are the terms of 

(delG p + G v). 

3) The third section contains the pressure adjacency terms (pi - pi+i). These 
values should be close to zero. 

4) The fourth section contains the velocity adjacency terms (vi - Vi+i). These 
values should be close to zero. 

The following are the primary criteria that impact the successful use of this 
technique for far field extrapolation. 

1) The complex response data (section two) of the two matrices, Ma-new and 
Ma, should agree with one another as closely as possible. 

2) The (-Ap + Bv) term (section one) of the Ma-new matrix should be 
minimized as much as possible since this represents a known physical parameter. 

3) The pressure adjacency and velocity adjacency terms (sections three and 
four) are expected to be small terms, but this condition does not require a heavy 
enforcement. Some differences in pressure and velocity values on adjacent 
patches are expected. This criteria was originally added to lessen the chance of 
skewing the pressure and velocity values to physically unrealizable states. 
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4) Since “zeros” are not likely to be obtained for seetions 1,3 or 4, it is 
important to eonsider the differenee (in dB) between the magnitudes in these 
seetions and the magnitudes of the eomplex response data at each element of the 
reeeive array (seetion 2). This differenee should be maximized as mueh as 
possible. 

The data presented below shows the two matriees, Ma-new and Ma, for a given set¬ 
up. This data presented has been optimized to a reasonable level (as examined later) and 
would serve to qualify the solved-for pressure and veloeities for the last step of this 
methodology, the step when the far field eomputations are aetually made. 

The four distinct sections are most easily observed in the Ma-new figure. The 
seetions 1 thru 4, as defined above, appear in suecession from the top of the figure to the 
bottom. The figure shows magnitudes expressed in dB; the colorbars to the right of the 
figures designate the eorrespondence between eolor and magnitude. The reddish band of 
colors correspond to the complex response data (seetion 2) and, as expeeted, are the 
highest magnitudes in the matrix. 



Figure 4 - Original and “pv solution” Augmented Field Point Matrices (Ma and Ma-new) 
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I. 


COMPUTE FIELD POINT RESPONSES 


This last step uses the pressure and veloeity data obtained and evaluated from the 
previous two steps to compute complex response at defined field points. This 
computation proceeds via the use of the Greens function integral and its gradient as has 
been accomplished in earlier steps when defining the augmented matrices. This familiar 
relationship is: 

p = (delG p + G v) 

As before, the p and v values represent the pressures and velocities at each patch 
location on the enclosing sphere. The G and delG terms represent the Greens function 
integral and its gradient obtained at specified field points with respect to the enclosing 
sphere. The selection of field points will typically be made to encircle the acoustic 
device on a designated plane allowing what are traditionally called XY, XZ, or YZ 
directivities. 

It is significant to note here that there is no distinction in this step whether defined 
field points are in the near field or in the far field. This technique takes a 
pressure/velocity defined enclosing sphere and combines its elemental traits via known 
physical relationships (Greens function integral and its gradient) to achieve pressure 
values at defined locations. Normalizing directivity data obtained via this technique for 
an increasing range will show the expected convergence to a single far field pattern of 
defined character. 
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V. SOLUTION SET USING METHODOLOGY 


This section takes the theory and methodology detailed in the previous sections to 
solve a real world problem. 

A, DEFINITION OF SET-UP 

In this case (and as supported by the referenced figures, Table 1 and the 
MATLAB coding provided), a horizontal transmitting line array of length 1.885 meters is 
operated at a frequency of approximately 4775 Hz. The expected near field/far field 
boundary for this device occurs at a distance of approximately 11.32 meters. This 
distance exceeds the physical size of the Acoustic Pressure Tank Facility (APTF), 11.1 
meters in length, where this device’s performance under temperature and pressure 
conditions are desired to be characterized. It is an ideal candidate to investigate the 
methodology detailed within this thesis. 

Figures 8 and 9 show the expected theoretical near field and far field responses, 
respectively, of the defined device. Figure 8 shows near field performance at a range of 
approximately 5.341 meters. This test distance is within the available range of separation 
distances that can be achieved within the APTF. The theoretical far field response was 
computed at a test distance of about 1131 meters. This places device performance 
characterization well into far field conditions. (Please note that the designation of 
element #6 on the near field data , Figure 8, corresponds to response in the XY plane. As 
will be described later, element #6 is the mid-point element of the 11 element vertical line 
array used in this effort.) 

A vertical near field array of approximate length 1.571 meters was used for this 
set-up. This near field array is made up of 11 equally spaced elements and was placed at 
a distance of 5.341 meters from the origin. (The horizontal transmitting array will be 
rotated about this same origin in the XY plane for this set-up. This origin location 
corresponds to the mid-point of the transmitting array.) The lengths specified for both 
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the near field array and for the transmitting array fit within the eonfines of the 3.81 
diameter APTF test vessel. 

The seleetion of 11 elements for the near field array was arrived a bit by 
experimentation and by a eoneerted effort to reduee the numbers of aequisition ehannels 
that ultimately would be required for this implementation (ie - 16 ehannel aequisition 
eards are fairly easy to proeure as off-the-shelf eomponents and would be suited to a 
future hardware implementation for this teehnique). 

Figure 7 shows a seetion of the enelosing CHIEF sphere that was used. The 
enclosing sphere was configured to have a radius of about 1.384 meters. The sphere 
section displayed in Figure 7 is made up of 84 individual patches. A symmetry order of 
180 was selected; this yields a total of 15,120 patches on the enclosing sphere surface. It 
is a part of this methodology to solve for a single pressure and a single velocity at each 
one of these patches. Thus, 30,240 “unknowns” must be solved for in this case. (As will 
be discussed later, the number of individual patches required on the enclosing surface is 
the primary limiting factor for this methodology. High frequency devices with large 
apertures will likely require an enclosing sphere of relatively large size yielding an 
extreme number of patches. This likely prohibits solution via this implementation on 
standard personal computing hardware and available COTS software.) 

B, INITIAL TEST RESULTS 

Figures 10 through 14, inclusive, show outputs obtained via the second step of the 
MATLAB coding. Figure 10 shows a comparison in absolute magnitude of the original 
data set (obtained via simulation in this example but would correspond to actual 
measured data in a real world implementation) and the hydrophone near field response 
obtained via the technique. These latter hydrophone responses were obtained using the 
solved for pressure and velocities from the technique. Figure 10 actually shows the 
difference between these two data sets. All differences are less than |0.015| dB showing 
an excellent match. 
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The remaining figures, 11 thru 14, allow the operator to assess the other 
eonditions of the methodology. Specifically, the values computed for the Helmholtz 
Integral relations and those corresponding to the pressure and velocity adjacency 
conditions need to be assessed as well to judge acceptance of the solved for pressures and 
velocities. In an ideal solution, these other conditions would be equal to zero. Since the 
least squares solution technique was utilized, a best fit to the defining equations was 
determined based upon near field array element data and assumed zero conditions for the 
Helmholtz Integral relations and pressure/velocity adjacency conditions. Is it not a true 
requirement to have these latter two conditions exactly equal to zero. It is, however, 
desirable to have these other conditions much smaller than the computed hydrophone 
responses. It is for this reason that dB levels are reported in Figures 11 through 14. 

Figure 11 shows a plot of the minimum hydrophone responses, maximum ZeroO 
and maximum Zerol conditions on a per angle basis. The nomenclature ZeroO and Zerol 
correspond to the “expected zeros” for the Helmholtz Integral relations and 
pressure/velocity adjacency conditions, respectively. Since the need here is to assess the 
solved for pressures and velocities, it is desirable to look at worst case conditions in these 
data sets. That is, the selection of minimum hydrophone response and maximum ZeroO 
and Zerol responses provides the smallest difference between the two data set conditions. 
(Reference Figure 3 - Original and “pv solution” Augmented Field Point Matrices (Ma 
andMa-new) for visualization of the ZeroO and Zerol conditions. ZeroO values are shown 
above the reddish band, the complex response data, while Zerol values are shown 
below.) 


Figures 12 and 13 augment Figure 11 by providing the actual computed 
differences described above. In this data set, a minimum difference on the order of about 
44 dB and 54 dB are observed for the ZeroO and Zerol conditions. These magnitudes, 
from experimentation, are likely to indicate that the solved for pressures and velocities 
are acceptable to use for far field extrapolation. The adjustment of the values for the 
alphaO and betaO terms will be looked at later to try and better this solution. 
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Figure 14 is perhaps the most useful of all of the figures in the Figure 11 through 
14 set. Figure 14 provides a “snapshof’ of all of the parameters discussed above 
(Helmholtz Integral Relations, Near Field Hydrophone Response, Pressure Adjacency 
and Velocity Adjacency conditions). It provides a means to directly compare these 
parameters with their originating values. It also allows the operator to make direct 
comparison amongst the Ma-new constituents (which are the parameters mentioned above). 
The color in this figure is representative of magnitude as shown in the colorbar scale 
provided. Thus, it is possible to roughly deduce at a glance, all of the differences as noted 
in Figure 11 through 13. 

Figures 15 through 18 used the solved for pressures and velocities obtained for 
alphaO=l and betaO=l (with the performance characteristics as identified in the previous 
figures) to solve for pressure at defined field points. These defined points correspond to 
the original near field points and points in the far field (which is the motivation for this 
technique). Figures 15 and 16 show that the solved for pressures and velocities yield near 
field pressures that have an extremely high degree of agreement with the original near 
field data set. Absolute differences of less than .008 dB were observed. 

Figures 17 and 18 show far field response as originally simulated and as 
extrapolated by this technique. The majority of the differences in these two data sets 
(Figure 18) are on the order of a few tenths of a dB with just a few outliers (these are 
addressed later). Thus, the technique has done an excellent job of matching expected far 
field performance with that obtained via the use of solved for pressures and velocities on 
an enclosing sphere. The technique is therefore viable. 

As mentioned above, the majority of the differences in the two data sets (Figure 
18) are on the order of a few tenths of a dB. There are, however, some difference values 
as high as about 4.3 dB. Closer examination of the locations of these large differences 
shows a correspondence with a null in the directional response. Thus, the technique did 
an excellent job of producing a far field response but missed representing to a high 
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accuracy the depth of the nulls in these far field responses. This will be looked at again 
once new values of alphaO and betaO are utilized. 

The values for alphaO and betaO are now changed to evaluate their effect on the 
solved for pressures and velocities and eventually their effeet on extrapolated far field 
performanee. 

C. SENSITIVITY OF SOLUTION TO VARIATION IN BETAO 

Figures 19 through 23 used values of alphaO=l and beta0=10. (Reeall that the 
betaO term emphasizes near field element data.) Figure 23 shows some improvement in 
the difference in level between the hydrophone response measured and its other zeroO 
and zerol eonstituents. These minimum differences are now about 74 and 64 dB, 
respectively. Pressure eomputations at the various field points were not eomputed. 

Figures 24 through 28 used values of alphaO=l and beta0=100. An improvement 
in difference levels was noted (97 dB and 84 dB), however, pressure computations were 
not undertaken. 

Figures 29 through 33 used values of alphaO=l and beta0=500. Differences of 
110 dB and 98 dB between hydrophone response and ZeroO and Zerol conditions, 
respectively, were noted. Figure 33 shows this same level of differenee by inspeetion. 

Despite ehanges in magnitude between these various Ma-new figures, an 
examination of Figure 33 and Figure 14 shows how the betaO term is foreing a larger 
difference between the hydrophone response data and the Helmholtz Integral relations 
and pressure/veloeity adjacency conditions. Each of these Ma-new figures has a dynamic 
range of 150 dB. Figure 33 shows how the hydrophone response data is being 
emphasized while the other parameters are being suppressed. 

Computation of pressure at the various field points was undertaken for the 
alpha0=l and beta0=500 conditions. Figures 34 through 37 show the resulting pressure 
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responses and simulated/teehnique derived differenees for these values. As before, the 
near field data eontinues to show solid agreement between simulated and teehnique 
derived values (Figures 34 and 35). Figures 36 and 37 show far field performanee traits. 
The noteworthy part of this data set is that the differenees between the two far field data 
sets are now less than 1 dB (down from 4.3 dB in Figure 18). Thus, this betaO value 
ehange has resulted in pressures and veloeities that provide a better far field fit. This is 
espeeially seen now in the better agreement afforded at the nulls of the far field response. 

A final betaO ehange to betaO = 1000 was undertaken. Figures 38 through 43 
show data obtained with alpha0=l and beta0=1000. As ean be seen on the far field 
differenee reported in Figure 43, there was little ehange in far field response provided by 
this. 


While far field agreement in this ease eould have essentially stopped with 
alpha0=l and beta0=l, the other values of betaO were examined to evaluate the 
teehnique. Inereases in betaO did lead to a better fit in the nulls of the far field, but this is 
likely of limited praetieal value. 

D, SENSITIVITY OF SOLUTION TO VARIATION IN ALPHAO 

In likewise fashion, to further evaluate the teehnique, ehanges in the value of 
alphaO (emphasizing the theoretieal portion of the teehnique) will now be undertaken. 

Figures 44 through 48 show data eorresponding to alpha0=10 and beta0=l. 

Figures 49 through 53 show data eorresponding to alpha0=100 and beta0=l. 

Figures 54 through 58 show data eorresponding to alpha0=100 and beta0=l. 

An examination of the Ma-new figures for eaeh of these eonditions shows that there 
is little to no ehange in level differenees between hydrophone response and 


34 



pressure/velocity adjacency conditions. The Helmholtz Integral relations are showing a 
suppression as alphaO increases. 

Figures 59 through 62 show pressure responses computed for the near field and 
far field cases. An examination of the far field difference figure (Figure 62) shows that 
very little was gained by increasing the alphaO term. The differences in the nulls that 
were noted for the alphaO=l and betaO=l far field computations still exist and are at the 
same approximate level. 

E, SUMMARY OF SENSITIVITY OF SOLUTION TO ALPHAO AND BETAO 

Table 1 summarizes the conditions examined in the previous examples. The 
following definition of terms applies to the data presented. 

alphaO: 

This is the weighting term to emphasize the theoretical Helmholtz Integral 
Relations. 

betaO: 


This is the weighting term to emphasize empirical data 

MAX |Near Field Hydro Diflj: 

This value represents the maximum difference between the measured hydrophone 
responses of all near field receive array elements and the values computed by the second 
step of the methodology. A low value represents good agreement and, therefore, a 
possible candidate for far field extrapolation (using the pressure and velocity values so 
obtained). 
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MIN min hydro minus max zeroO: 

This value represents the minimum differenee between the minimum near field 
hydrophone data and the maximum zeroO data (theoretieal Helmholtz relations) obtained 
at eaeh angular position. A large magnitude here also signifies that the eonditions are 
favorable for far field extrapolation. 

MAX min hydro minus max zeroO: 

This value represents the maximum differenee between the minimum near field 
hydrophone data and the maximum zeroO data (theoretieal Helmholtz relations) obtained 
at eaeh angular position. This is largely for informational purposes only. 

MIN min hydro minus max zerol: 

This value represents the minimum differenee between the minimum near field 
hydrophone data and the maximum zerol data (pressure and veloeity equivaleney for 
adjaeent patehes) obtained at eaeh angular position. A large magnitude here also 
signifies that the eonditions are favorable for far field extrapolation. 

MAX min hydro minus max zerol: 

This value represents the maximum differenee between the minimum near field 
hydrophone data and the maximum zerol data (theoretieal Helmholtz relations) obtained 
at eaeh angular position. This is also largely for informational purposes only. 

MAX |Far Field Difference!: 

This value represents the maximum differenee between the methodology- 
eomputed far field response and a measured (or in this ease, simulated) far field response 
obtained at eaeh angular position. A low value indieates that the methodology produeed 
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a reasonable approximation to aetual far field performanee. A high value may indicate an 
overall poor match or perhaps just a poor representation at certain angles (most likely in 
the nulls). It is advantageous to look at the overall pattern obtained to further define the 
severity of a high value. (Please note, however, that all constituents to obtain this data 
value would not typically be available in practice (ie - the purpose of this thesis was to 
achieve a far field formulation when only near field test conditions were available). 
These results were available in this thesis effort and therefore have been examined.) 

The initial equal weighting of the alphaO and betaO terms (by setting both equal to 
1) for this example gave excellent results with little gained by any other adjustments 
(other than a bit better definition of the response in far field nulls). This is not necessarily 
likely to be accurate in all cases. Where there is some uncertainty in the measured data 
acquired, a better fit may be obtained by emphasizing the theoretical portion of the 
technique. This will be examined in the next section when error in range and positioning 
of the near field array are examined. 
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MAX 


MIN 


MAX 


alphaO betaO 


INear Field Hydro Diff| 


min_hydro_minus_max_zeroO 


min_hydro_minus_max_zeroO 


(dB) 


(dB) 


(dB) 


1 

1 

0.010961319 

54 . 0 

1 

10 

0 . 001350528 

74.2 

1 

100 

0.000097393 

96.7 

1 

500 

0.000008493 

110.5 

1 

1000 

0.000009353 

116.5 

10 

1 

0.011399652 

71 . 7 

100 

1 

0 . 011405154 

91 . 7 

500 

1 

0 . 011405208 

105 . 7 


93.9 

113.9 
132 . 6 
146.2 
152 . 1 
111 . 4 

131 . 4 

145 . 4 


MIN MAX MAX 

alphaO betaO min_hydro_minus_max_zerol min_hydro_minus_max_zerol | Far Field Difference! 




(dB) 

(dB) 

(dB) 

1 

1 

43.8 

89.2 

4.30 

1 

10 

63.8 

107 . 6 

n/a 

1 

100 

84.2 

125.2 

n/a 

1 

500 

97 . 9 

138.9 

0.95 

1 

1000 

103.8 

144 . 8 

0.90 

10 

1 

43.6 

88.7 

n/a 

100 

1 

43.6 

88.6 

n/a 

500 

1 

43.6 

88.6 

4.00 



Table 1 - 

Solution Set Methodology Data Summary 
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VI. ERROR ANALYSIS - EFFECT OF POSITIONAL UNCERTAINTY 


A, INTRODUCTION 

This section addresses positional uncertainties in the placement of the near field 
array and its resulting effect on the extrapolated far field response using the described 
methodology. The examples provided herein were computed using the same device 
configurations, acoustic parameters and CHIEF set-ups detailed in the previous Solution 
Set using Methodology section. In this section, however, the alphaO and betaO terms 
were established at unity each and were not allowed to vary. 

Positional uncertainty was defined to correspond to an inexact determination of 
range, vertical placement or angular orientation of the near field array with respect to the 
transmitting array’s origin. The data for this section was obtained by varying one of the 
aforementioned parameters and then running through the computation steps in the 
methodology to achieve extrapolated far field response performance. The assessment of 
the error so induced is primarily based upon the examination of the far field responses 
and their differences. Near field responses (showing excellent agreement throughout the 
various error conditions) and figures of the Ma-new matrices are provided for 
completeness in this section. 

B, RANGE ERROR 

Figures 63 through 86, inclusive, examine the effect of intentionally induced 
errors in range. For these figures, the simulated data used range values with error while 
the thesis technique used the nominal range magnitudes (ie: range without error). 
Various range errors were utilized to assist in evaluating the sensitivity of the technique 
to this parameter. These are summarized on the following page. 
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Figures 63 through 68 


-2*lambda (-0.628 meters) 


Figures 69 through 74 +2*lambda (+0.628 meters) 

Figures 75 through 80 -1 em 

Figures 81 through 86 +1 em 

For referenee, Figures 10 through 18 eorrespond to a range error equal to zero 
(and with alpha0=l and beta0=l). 

As expeeted, the larger range errors provided the most error in the extrapolated far 
field response as shown in the far field differenee data reported in Figures 68 and 74. 
Even in these extreme eases (+/- 0.628 meters on a defined range of about 5.354 meters), 
however, the far field response (negleeting the response in the nulls of the response) 
differed from the zero range error ease on average no more than about 0.5 dB. These 
results were a bit unexpeeted. 

The range errors of +/- 1 em were eomputed in this data set to better represent 
range errors that are more likely to be encountered within the test facility (APTF) . 
These error conditions (Figures 74 and 80) show very little difference from the zero range 
error case. Thus, expected uncertainties in the range (to +/- 1 cm and for the set-up 
conditions examined) are not deemed to be a significant source of error in extrapolating 
far field response via the technique presented. 

C. VERTICAL PLACEMENT ERROR 

Figures 87 through 104, inclusive, examine the effect of intentionally induced 
errors in the vertical placement of the receive array. In like fashion as for the range error 
case, the simulated data used the error condition while the thesis technique used the 
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nominal values (ie; without error). Various errors were utilized to assist in evaluating the 
sensitivity of the technique to this parameter. These are summarized below; 

Figures 87 through 92 +lambda/2 (+0.157 meters) 

Figures 93 through 98 +lambda (+0.314 meters) 

Figures 99 through 104 -lambda/2 (-0.157 meters) 

The selection of these error values as multiples of lambda/2 stems from the 
configuration of the near field receive array itself. In the example that has been 
discussed, the spacing of the elements in the near field array is lambda/2. The worst error 
condition that was deemed likely to occur in practice was considered to be an offset in the 
placement of the array by one element. A 2 element offset (+lambda in this example) 
was selected to further assess the technique. This latter vertical error is, however, not 
likely to occur in practice for the inter-element spacings specified in this example. As for 
the range error case. Figures 10 through 18 can be referenced for comparison to a zero 
error case. 

Again, as expected, the larger vertical placement errors provided the most error in 
the extrapolated far field response as shown in the difference data reported in Figure 98. 
For the more typical error cases (vertical errors of +/- lambda/2), the differences in far 
field response were comparable with those obtained for the zero error case. None of the 
vertical placement errors accounted for more than perhaps a 0.25 dB error in the final 
reported far field response (again, neglecting the response in the nulls of the response). 

Thus, uncertainties in the vertical placement of the near field array (to within an 
element spacing and for the set-up conditions examined) are not deemed to be a 
significant source of error in extrapolating far field response via the technique presented. 
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D, ANGULAR ORIENTATION ERROR 


Figures 105 through 110, inclusive, examine the effect of intentionally induced 
errors in the angular orientation (offset) of the receive array with respect to the 
transmitting array. In like fashion as for the previous two cases, the simulated data used 
the error condition while the thesis technique used the nominal values (ie: without error). 
A single value of error equivalent to a 10 degree offset was examined in this data set. 

Figure 110 shows the difference between the simulated data set and that obtained 
from the thesis technique. Upon initial inspection, this figure seems to indicate gross 
errors caused by the 10 degree displacement. Differences of close to 30 dB are noted in 
this figure. Upon closer examination of the defining constituents of the difference figure 
(the data presented in Figurel09), it is readily seen that the shape and structure of the two 
far field responses are essentially the same. The introduction of the 10 degree angular 
error has shifted the response accordingly and caused the large difference errors. While 
this results in large difference errors (Figure 110), this condition is easily post-processed 
and does not represent any substantial source of error or a deterrent to the use of the 
thesis methodology. 

Thus, as for the other cases, uncertainties in the angular orientation of the near 
field array are not deemed to be a significant source of error in extrapolating far field 
response via the technique presented. This statement is likely to hold true for all set-ups 
and not just for the unique parameters looked at in this example set. 

E, ERROR SUMMARY COMMENTS 

Table 2 summarizes the error conditions examined in the previous examples in 
this section. Please reference the previous section. Solution Set using Methodology, for 
definitions of the nomenclatures utilized. 
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For the given example set-up, errors in range, vertieal plaeement and angular 
orientation of the near field reeeive array (when taken singularly and not taken in any 
kind of an extreme sense) indueed very little error in the extrapolated far field response 
(again, negleeting the nulls of the response) Attention to error quantifieation in terms of 
some lambda value (or fraetional value) was given in an effort to assist in extrapolating 
this performanee data to other set-up eonditions. While this ability to extrapolate 
performanee traits is likely aehievable, the ability to position deviees within the same 
lambda toleranee will be hampered as operating frequeneies inerease (and wavelengths 
eorrespondingly deerease). 

The overall eonelusion provided, however, is that the teehnique worked extremely 
well at eharaeterizing far field response even when physieal set-up parameters were 
known to be in error. 
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alphaO=betaO=l 


MAX 


MIN 


MAX 


error condition 
no error 

range error: -2*lambda 
range error: +2*lambda 
range error: -1 cm 

range error: +1 cm 
vert rec array error: 

+lambda/2 

vert rec array error: ilambda 
vert rec array error: - 
lambda/2 

azimuthal error: +10 degrees 


error condition 
no error 

range error: -2*lambda 
range error: +2*lambda 
range error: -1 cm 

range error: +1 cm 
vert rec array error: 

+lambda/2 

vert rec array error: llambda 
vert rec array error: - 
lambda/2 

azimuthal error: +10 degrees 


I Near Field Hydro Diff| 
(dB) 

0.010961319 
0 . 015705251 
0 . 010182845 
0 . 011011699 
0 . 010911516 

0 . 019242789 
0 . 023584640 

0 . 019242789 
0 . 011221058 

MIN 

min_hydro_minus_max_zerol 

(dB) 

43.80 
40.43 

41 . 85 
43.83 
43 . 78 

35.09 

29.10 

35.09 

43.80 


min_hydro_minus_max_zeroO 

(dB) 

53.95 

49.88 
53.68 

53.96 
53.95 

50.60 

40 . 73 

50.60 
53 . 95 

MAX 

min_hydro_minus_max_zerol 

(dB) 

89.16 

83.27 

92 . 73 
89.05 
89.26 

87.88 
81.83 

87.88 

89.16 


min_hydro_minus_max_zeroO 

(dB) 

93.86 

88.29 

93.62 

93.78 

93.93 

89.97 

85.23 

89.97 

93.86 

MAX 

I Far Field Difference! 
(dB) 

4.30 

12.60 
5.20 
4.40 
4.25 

3.00 

2.60 

3.00 

30.00 


Table 2 - Effect of Positional Uncertainty Data Summary 
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VII. ADDITIONAL SOLUTION SETS 
A, INTRODUCTION 

This section addresses the effectiveness of the technique when the near field array 
is placed closer to the device under test. The same transmitting array, 1.885 meters in 
length and operating at a frequency of about 4775 Hz, is used in all cases. The cases 
presented in this section are summarized below: 



Near Field Array Location 

alphaO 

betaO 

Figures 111 to 113 

2.082 meters 

1 

1000 

Figures 114 to 116 

1.466 meters 

1 

1000 

Figures 117 to 119 

1.466 meters 

1 

10,000 

*Figures 120 to 122 

1.466 meters 

1 

1000 


* - Figures 120 to 122 represent data taken with a different configuration near field 
array. The array used for this data set was about 2.199 meters long and was composed of 
15 equally spaced elements. 

The data presented here was provided more for academic reasons than for any 
required need. The location of the near field array at a range of 5.354 meters in the 
previous sections adequately met the restrictions imposed by the dimensions of the APTF 
test vessel. Thus, the test set-ups specified in this section may not be required, but they 
do serve to address the capabilities of the technique. 

B. NEAR-FIELD ARRAY RANGE OF 2,082 METERS 

Figures 111 through 113 used a near field array range of 2.082 meters. The near 
field data sets (simulated and from the technique) matched very well; the far field data 
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sets also matched extremely well. It is interesting to note that the near field data actually 
now shows a local null at the same angular location where the far field data actually 
peaks. 


Please note that the enclosing sphere required for these data sets was also moved 
to a closer proximity to the device under test. Since the radius of this enclosing sphere 
was decreased and yet the frequency of operation stayed the same, fewer numbers of 
patches were required in the vertical direction and a lower symmetry order sufficed for 
the data set. The lower symmetry order (since it is directly related to the number of 
discrete angles in the responses) accounts for the choppiness of these responses when 
compared to the previous solution sets’ figures. 

C. NEAR-FIELD ARRAY RANGE OF 1,466 METERS 

Figures 114 to 122 all used a near field array range of 1.466 meters. Figures 114 
through 116 and Figures 117 through 119 used the same alphaO values (one) but differing 
betaO values of 1000 and 10,000, respectively. The far field response of Figure 115 
showed a very reasonable fit to the simulated data along the main beam of the device, but 
did not have as good a match in the side lobes. The data for Figure 118 was computed 
with a higher betaO value to attempt to correct this side lobe problem. The Ma-new data for 
these two conditions (Figures 116 and 119) would lead one to believe that the higher 
betaO value would give a more accurate description of the devices actual far field 
response. This was not found to be the case. The far field response for these two 
conditions (Figures 115 and 118) looked to be very comparable. 

After some experimentation, it was observed that increasing the length of the near 
field array yielded an improvement in side lobe performance of the device under test. 
Figures 120 to 122 show response with a near field array that has been lengthened by 
2*lambda to 2.199 meters (originally this array was about 1.571 meters in length). The 
far field response (Figure 121) shows an improved side lobe performance over the 
previous data sets using the same near field array range of 1.466 meters. 
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Please note that the number of elements in the 2.199 meter near field array was 
inereased from 11 to 15 in an effort to maintain lambda/2 spaeing. The improvement in 
far field response observed, however, was due to its length inerease and not due to the 
inerease in the number of its eomprising elements. This latter characteristic was arrived 
at empirically but is not supported with any data in this write-up. 

D, SUMMARY 

Based upon the data sets presented in this section, it is desirable to quantify a 
minimum near field array length and a minimum test distance for a given test set-up to 
provide some assurance that proper far field performance is extrapolated. Using 
empirical results, it is believed that a near field array with a length that exceeds the 
maximum dimension of the device under test placed at a separation distance comparable 
to its length could achieve this condition. This latter hypothesis is speculated but is not 
stated as proof in this thesis. It is however, supported by the data presented in Figures 
123 through 125 (with the transmitting array, near field array and range all specified at 
1.885meters). 
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VIII. SUMMARY 


A, INTRODUCTION 

A technique to extrapolate far field performance of an acoustic device from a data 
set obtained in the near field was demonstrated. The technique successfully married the 
Helmholtz integral relationships (A and B matrices and pressure computations utilizing 
the integration of the Green’s function and its gradient) and expected physical parameters 
with near field acquired data to yield an over determined system. Solution via a least 
squares technique provided pressure and normal velocity terms for each patch of the 
enclosing sphere surface. This characterization of the enclosing sphere’s surface allowed 
extrapolation to any range (near or far field) when Greens’ function and its gradient 
integration values appropriate to that range are utilized. 

B, SOLUTION SETS 

An example solution set was provided that demonstrated the ability of the 
technique to produce extremely accurate results. Figures 4 and 5 below show the data 
computed for the near field and far field, respectively, for the defined acoustic device (a 
1.885 meter long linear array operating at about 4775 Hz). (Figures 4 and 5 represent the 
same data set shown in Figures 34 and 36.) 
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NEAR FIELD DATA (Rec Ele #6) with k = 20; range = 5,3541 meters; array length = 1,885 meters 



Figure 5 - Near Field Data - computed and simulated 


FAR FIELD DATA (on axis) with k = 20; range = 1130,9734 meters; array length = 1,885 meters 



Figure 6 - Far Field Data - computed and simulated 


The example set-up was specifically designed for a minimal set of elements in the 
near field array (eleven). This minimal set was utilized in an attempt to ease the 
adaptation of the technique to a future hardware implementation (ie: 16 channel 
acquisition cards of sufficient dynamic range and frequency of operation are readily 
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available now; the eleven near field array would require only a single eard for 
implementation.) 

C. SUMMARY 

It was desired to verily the teehnique with aetual measurement data, however, the 
present eonliguration of the Aeoustie Pressure Ta nk Facility’s measurement system ( 4 
channel ; VXI-based) precluded a reasonably fast or convenient way to measure the 
complex response of an eleven element near field array. This latter array also does not 
presently exist, but likely could have been conjured up within existing USRD inventories 
if so required. It is anticipated that this technique will be looked at with an updated 
acquisition system and perhaps a specially designed array at some future date. 

An error study was undertaken with the simulated data sets. This error study 
examined errors in range, vertical placement and angular orientation of the near field 
array. This study confirmed that the technique (at least for the given set-up) was 
extremely robust and did not require exacting tolerances to achieve qualifiable results. 
This was especially gratifying for the test environment in which the technique is 
ultimately intended to be used. 

The technique will not, however, be appropriate for all devices and all operating 
frequencies. The reliance on COTS computing hardware and commercial software 
products does limit the size of the matrices that can be reasonably handled. As 
mentioned previously, the solution set example used a relatively low frequency of about 
4775 Hz but yet still required the solution of 30,240 “unknowns” for its 15,120 individual 
patches. 

In summary, the technique does present itself favorably and is highly likely to be 
looked at for future implementation at USRD. 
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FAR FIELD EXTRAPOLATION TECHNIQUE USING 
CHIEF ENCLOSING SPHERE DEDUCED PRESSURES AND VELOCITIES 

FIGURES 
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Figure 7 - Enclosing Sphere Section (numv =84; symmetry=180) 


NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 8 - Theoretical Near Field Response of Finear Array 
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FAR FIELD DATA with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 9 - Theoretical Far Field Response of Linear Array 
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Figure 10 - Near Field Flydrophone Difference Data 
alpha0=l; betaO = 1 


Ma-new - min(HYDROPHONE), max(Zero0), max(Zero1); alphaO = 1 ; alpha = 1.9789 ; betaO = 1 ; beta = 998.5472 



Figure 11 - Ma-new Components 
alphaO=l; betaO = 1 
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Ma-new [ min(HYDROPHONE) - max(ZeroO) ] (dB); alphaO = 1 ; betaO = 1 



Figure 12 - Ma-new [min(Hydrophone) - max(ZeroO)] Data 
alphaO=l; betaO = 1 


Ma-new [ min(HYDROPHONE) - max(Zero1) ] (dB); alphaO = 1 ; betaO = 1 



Figure 13 - Ma-new [min(Hydrophone) - max(Zerol)] Data 
alphaO=l; betaO = 1 
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Figure 14 - Ma-new and Ma (original) Data Sets 
alphaO=l; betaO = 1 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 15 - Near Field Data - computed and simulated 
alpha0=l; betaO = 1 


NF DATA Difference (Rec Ele #6) with k = 20 ; range = 5.3541 meters 



Figure 16 - Near Field Data Difference 
alphaO=l; betaO = 1 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 17 - Far Field Data - computed and simulated 
alpha0=l; betaO = 1 


FF DATA Difference(on axis) with k = 20; range = 1130.9734 meters 
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Figure 18 - Far Field Data Difference 
alpha0=l; betaO = 1 
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Figure 19 - Near Field Flydrophone Difference Data 
alphaO=l; betaO =10 


Ma-new- min(HYDROPHONE), max(ZeroO), max(Zero1); alphaO = 1 ; alpha = 1.9789 ; betaO = 10 ; beta = 9985.4718 



Figure 20 - Ma-new Components 
alpha0=l; betaO =10 
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Ma-new [ min(HYDROPHONE) - max(ZeroO) | (dB); alphaO = 1 ; betaO = 10 



Figure 21 - Ma-new [min(Hydrophone) - max(ZeroO)] Data 
alpha0=l; betaO =10 


Ma-new [ min(HYDROPHONE) - max(Zero1) ] (dB); alphaO = 1 ; betaO = 10 



Figure 22 - Ma-new [min(Hydrophone) - max(Zerol)] Data 
alpha0=l; betaO =10 


62 





















Figure 23 - Ma-new and Ma (original) Data Sets 
alphaO=l; betaO =10 
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HYDROPHONE DATA - dB[(compufed-actual)]; alphaO = 1 ; betaO = 100 
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Figure 24 - Near Field Flydrophone Difference Data 
alphaO=l; betaO = 100 


Ma-new- min(HYDROPHONE), max(ZeroO), max(Zero1); alphaO = 1 ; alpha = 1.9789 ; betaO = 100 ; beta = 99854.7185 



Figure 25 - Ma-new Components 
alphaO=l; betaO = 100 
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Figure 26 - Ma-new [min(Hydrophone) - max(ZeroO)] Data 
alphaO=l; betaO = 100 


Ma-new [ min(HYDROPHONE) - max(Zero1) ] (dB); alphaO = 1 ; betaO = 100 
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Figure 27 - Ma-new [min(Hydrophone) - max(Zerol)] Data 
alphaO=l; betaO = 100 
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Figure 28 - Ma-new and Ma (original) Data Sets 


alphaO=l; betaO = 100 
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HYDROPHONE DATA - dB[(compufecl-actual)]; alphaO = 1 ; betaO = 500 



Figure 29 - Near Field Flydrophone Difference Data 
alpha0=l; betaO = 500 


Ma-new- min(HYDROPHONE), max(ZeroO), max(Zero1); alphaO = 1 ; alpha = 1.9789 ; betaO = 500 ; beta = 499273.5925 



Figure 30 - Ma-new Components 
alpha0=l; betaO = 500 
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Figure 31 - Ma-new [min(Hydrophone) - max(ZeroO)] Data 
alphaO=l; betaO = 500 


Ma-new [ min(HYDROPHONE) - max(Zero1) ] (dB); alphaO = 1 ; betaO = 500 



Figure 32 - Ma-new [min(Hydrophone) - max(Zerol)] Data 
alpha0=l; betaO = 500 
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Figure 33 - Ma-new and Ma (original) Data Sets 


alphaO=l; betaO = 500 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 
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Figure 34 - Near Field Data - computed and simulated 
alpha0=l; betaO = 500 

8 

X 10 '® NF DATA Difference (Rec Ele #6) with k = 20 ; range = 5.3541 meters 
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Figure 35 - Near Field Data Difference 
alphaO=l; betaO = 500 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 36 - Far Field Data - computed and simulated 
alpha0=l; betaO = 500 



Figure 37 - Far Field Data Difference 
alphaO=l; betaO = 500 
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HYDROPHONE DATA - dB[(computed-actual)]; alphaO = 1 ; betaO = 1000 



Figure 38 - Near Field Flydrophone Difference Data 
alpha0=l; betaO = 1000 


Ma-new - from pv solution 


Ma - original 



50 100 150 


50 100 150 


Figure 39 - Ma-new and Ma (original) Data Sets 
alpha0=l; betaO = 1000 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 40 - Near Field Data - computed and simulated 
alpha0=l; betaO = 1000 


X 10 '® NF DATA Difference (Rec Ele #6) with k = 20 ; range = 5.3541 meters 
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Figure 41 - Near Field Data Difference 
alphaO=l; betaO = 1000 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 42 - Far Field Data - computed and simulated 
alpha0=l; betaO = 1000 


FF DATA Difference(on axis) with k = 20; range = 1130.9734 meters 



Figure 43 - Far Field Data Difference 
alphaO=l; betaO = 1000 
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HYDROPHONE DATA - dB[(computed-actual)]; alphaO = 10 ; betaO = 1 



Figure 44 - Near Field Flydrophone Difference Data 
alpha0=10; betaO = 1 


Ma-new- min(HYDROPHONE), max(ZeroO), max(Zero1); alphaO = 10 ; alpha = 19.7892 ; betaO = 1 ; beta = 998.5472 



Figure 45 - Ma-new Components 
alpha0=10; betaO = 1 
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Ma-new [ min(HYDROPHONE) - max(ZeroO) | (dB); alphaO = 10 ; betaO = 1 



Figure 46 - Ma-new [min(Hydrophone) - max(ZeroO)] Data 
alpha0=10; betaO = 1 


Ma-new [ min(HYDROPHONE) - max(Zero1) | (dB); alphaO = 10 ; betaO = 1 



Figure 47 - Ma-new [min(Hydrophone) - max(Zerol)] Data 
alpha0=10; betaO = 1 
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Figure 48 - Ma-new and Ma (original) Data Sets 
alpha0=10; betaO = 1 
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HYDROPHONE DATA - dB[(compufecl-actual)]; alphaO = 100 ; betaO = 1 



Figure 49 - Near Field Flydrophone Difference Data 
alpha0=100; betaO = 1 


Ma-new- min(HYDROPHONE), max(ZeroO), max(Zero1); alphaO = 100 ; alpha = 197.8919 ; betaO = 1 ; beta = 998.5472 



Figure 50 - Ma-new Components 
alpha0=100; betaO = 1 
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Ma-new [ min(HYDROPHONE) - max(ZeroO) ] (dB); alphaO = 100 ; betaO = 1 



Figure 51 - Ma-new [min(Hydrophone) - max(ZeroO)] Data 
alpha0=100; betaO = 1 


Ma-new [ min(HYDROPHONE) - max(Zero1) ] (dB); alphaO = 100 ; betaO = 1 



Figure 52 - Ma-new [min(Hydrophone) - max(Zerol)] Data 
alpha0=100; betaO = 1 
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Figure 53 - Ma-new and Ma (original) Data Sets 


alpha0=100; betaO = 1 
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Figure 54 - Near Field Flydrophone Difference Data 
alpha0=500; betaO = 1 


Ma-new- min(HYDROPHONE), max(ZeroO), max(Zero1); alphaO = 500 ; alpha = 989.4593 ; betaO = 1 ; beta = 998.5472 



Figure 55 - Ma-new Components 
alpha0=500; betaO = 1 
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Ma-new [ min(HYDROPHONE) - max(ZeroO) ] (dB); alphaO = 500 ; betaO = 1 
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Figure 56 - Ma-new [min(Hydrophone) - max(ZeroO)] Data 
alpha0=500; betaO = 1 


90 


Ma-new [ min(HYDROPHONE) - max(Zero1) ] (dB); alphaO = 500 ; betaO = 1 
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Figure 57 - Ma-new [min(Hydrophone) - max(Zerol)] Data 
alpha0=500; betaO = 1 
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Figure 58 - Ma-new and Ma (original) Data Sets 
alpha0=500; betaO = 1 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 59 - Near Field Data - computed and simulated 
alpha0=500; betaO = 1 
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Figure 60 - Near Field Data Difference 
alpha0=500; betaO = 1 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 61 - Far Field Data - computed and simulated 
alpha0=500; betaO = 1 


FF DATA Difference(on axis) with k = 20; range = 1130.9734 meters 
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Figure 62 - Far Field Data Difference 
alpha0=500; betaO = 1 
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HYDROPHONE DATA - dB[(computed-actual)]; alphaO= 1 ; betaO= 1 



Figure 63 - Near Field Flydrophone Difference Data 
alphaO=l; betaO=l; range error = -2 * lambda 



Figure 64 - Ma-new and Ma (original) Data Sets 
alphaO=l; betaO^l; range error = -2 * lambda 


86 




































































NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 65 - Near Field Data - computed and simulated 
alpha0=l; betaO^l; range error = -2 * lambda 
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Figure 66 - Near Field Data Difference 
alphaO=l; betaO=l; range error = -2 * lambda 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 67 - Far Field Data - computed and simulated 
alpha0=l; beta0=l; range error = -2 * lambda 


FF DATA Difference(on axis) with k = 20; range = 1130.9734 meters 



Figure 68 - Far Field Data Difference 
alpha0=l; beta0=l; range error = -2 * lambda 
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HYDROPHONE DATA - dB[(computed-actual)]; alphaO= 1 ; betaO= 1 



Figure 69 - Near Field Flydrophone Difference Data 
alphaO=l; betaO=l; range error = +2 * lambda 



Figure 70 - Ma-new and Ma (original) Data Sets 
alpha0=l; beta0=l; range error = +2 * lambda 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 71 - Near Field Data - computed and simulated 
alpha0=l; beta0=l; range error = +2 * lambda 


NF DATA Difference (Rec Ele #6) with k = 20 ; range = 5.3541 meters 



Figure 72 - Near Field Data Difference 
alpha0=l; beta0=l; range error = +2 * lambda 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 73 - Far Field Data - computed and simulated 
alpha0=l; beta0=l; range error = +2 * lambda 
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Figure 74 - Far Field Data Difference 
alpha0=l; beta0=l; range error = +2 * lambda 
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Figure 75 - Near Field Flydrophone Difference Data 
alpha0=l; beta0=l; range error = -1 cm 



Figure 76 - Ma-new and Ma (original) Data Sets 
alpha0=l; beta0=l; range error = -1 cm 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 77 - Near Field Data - computed and simulated 
alpha0=l; beta0=l; range error = -1 cm 


NF DATA Difference (Rec Ele #6) with k = 20 ; range = 5.3541 meters 



Figure 78 - Near Field Data Difference 
alpha0=l; beta0=l; range error = -1 cm 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 79 - Far Field Data - computed and simulated 
alpha0=l; beta0=l; range error = -1 cm 


FF DATA Difference(on axis) with k = 20; range = 1130.9734 meters 



Figure 80 - Far Field Data Difference 
alpha0=l; beta0=l; range error = -1 cm 
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Figure 81 - Near Field Flydrophone Difference Data 
alpha0=l; beta0=l; range error = +1 cm 



Figure 82 - Ma-new and Ma (original) Data Sets 


alpha0=l; beta0=l; range error = +1 cm 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 83 - Near Field Data - computed and simulated 
alpha0=l; beta0=l; range error = +1 cm 


NF DATA Difference (Rec Ele #6) with k = 20 ; range = 5.3541 meters 



Figure 84 - Near Field Data Difference 
alpha0=l; beta0=l; range error = +1 cm 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 85 - Far Field Data - computed and simulated 
alpha0=l; beta0=l; range error = +1 cm 
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Figure 86 - Far Field Data Difference 
alpha0=l; beta0=l; range error = +1 cm 
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HYDROPHONE DATA - dB[(computed-actual)]; alphaO= 1 ; betaO= 1 



Figure 87 - Near Field Flydrophone Difference Data 
alphaO=l; betaO=l; vertical receive array error = +lambda/2 



Figure 88 - Ma-new and Ma (original) Data Sets 
alphaO=l; betaO=l; vertical receive array error = +lambda/2 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 89 - Near Field Data - computed and simulated 
alpha0=l; beta0=l; vertical receive array error = +lambda/2 


X 10 '^ NF DATA Difference (Rec Ele #6) with k = 20 ; range = 5.3541 meters 
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Figure 90 - Near Field Data Difference 
alpha0=l; beta0=l; vertical receive array error = +lambda/2 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 91 - Far Field Data - computed and simulated 
alpha0=l; beta0=l; vertical receive array error = +lambda/2 


FF DATA Difference(on axis) with k = 20; range = 1130.9734 meters 



Figure 92 - Far Field Data Difference 
alpha0=l; beta0=l; vertical receive array error = +lambda/2 
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Figure 93 - Near Field Flydrophone Difference Data 
alphaO=l; betaO=l; vertical receive array error = +lambda 



Figure 94 - Ma-new and Ma (original) Data Sets 


alphaO=l; betaO=l; vertical receive array error = +lambda 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 95 - Near Field Data - computed and simulated 
alpha0=l; beta0=l; vertical receive array error = +lambda 


NF DATA Difference (Rec Ele #6) with k = 20 ; range = 5.3541 meters 



Figure 96 - Near Field Data Difference 
alpha0=l; beta0=l; vertical receive array error = +lambda 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 97 - Far Field Data - computed and simulated 
alpha0=l; beta0=l; vertical receive array error = +lambda 


FF DATA Difference(on axis) with k = 20; range = 1130.9734 meters 
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Figure 98 - Far Field Data Difference 
alpha0=l; beta0=l; vertical receive array error = +lambda 
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HYDROPHONE DATA - dB[(computed-actual)]; alphaO= 1 ; betaO= 1 



Figure 99 - Near Field Flydrophone Difference Data 
alphaO=l; betaO=l; vertical receive array error = -lambda/2 



Figure 100 - Ma-newand Ma (original) Data Sets 
alpha0=l; beta0=l; vertical receive array error = -lambda/2 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 101 - Near Field Data - computed and simulated 
alpha0=l; beta0=l; vertical receive array error = -lambda/2 


X 10 '^ NF DATA Difference (Rec Ele #6) with k = 20 ; range = 5.3541 meters 
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Figure 102 - Near Field Data Difference 
alpha0=l; beta0=l; vertical receive array error = -lambda/2 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 103 - Far Field Data - computed and simulated 
alpha0=l; beta0=l; vertical receive array error = -lambda/2 


FF DATA Difference(on axis) with k = 20; range = 1130.9734 meters 



Figure 104 - Far Field Data Difference 
alpha0=l; beta0=l; vertical receive array error = -lambda/2 
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HYDROPHONE DATA - dB[(computed-actual)]; alphaO= 1 ; betaO= 1 



Figure 105 - Near Field Flydrophone Difference Data 
alpha0=l; beta0=l; azimuthal error = +10 degrees 



Figure 106 - Ma-newand Ma (original) Data Sets 
alpha0=l; beta0=l; azimuthal error = +10 degrees 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 5.3541 meters; array length = 1.885 meters 



Figure 107 - Near Field Data - computed and simulated 
alpha0=l; beta0=l; azimuthal error = +10 degrees 


NF DATA Difference (Rec Ele #6) with k = 20 ; range = 5.3541 meters 
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Figure 108 - Near Field Data Difference 
alpha0=l; beta0=l; azimuthal error = +10 degrees 
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FAR FIELD DATA (on axis) with k = 20: range = 1130.9734 meters; array length = 1.885 meters 



Figure 109 - Far Field Data - computed and simulated 
alpha0=l; beta0=l; azimuthal error = +10 degrees 


FF DATA Difference(on axis) with k = 20; range = 1130.9734 meters 



Figure 110 - Far Field Data Difference 
alpha0=l; beta0=l; azimuthal error = +10 degrees 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 2.0817 meters; array length = 1.885 meters 



Figure 111 - Near Field Data - computed and simulated 
Using 2.082 meter range, 1.571 meter long, 11 element near field array data 


FAR FIELD DATA (on axis) with k = 20; range = 1130.9734 meters; array length = 1.885 meters 



Figure 112 - Far Field Data - computed and simulated 
Using 2.082 meter range, 1.571 meter long, 11 element near field array data 
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Figure 113 - Ma-new and Ma (original) Data Sets 
Using 2.082 meter range, 1.571 meter long, 11 element near field array data 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 14658 meters; array length = 1.885 meters 



Figure 114 - Near Field Data - computed and simulated 
Using 1.466 meter range, 1.571 meter long, 11 element near field array data 


FAR FIELD DATA (on axis) with k = 20; range = 1130.9734 meters; array length = 1.885 meters 



Figure 115 - Far Field Data - computed and simulated 
Using 1.466 meter range, 1.571 meter long, 11 element near field array data 
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Figure 116 - Ma-new and Ma (original) Data Sets 
Using 1.466 meter range, 1.571 meter long, 11 element near field array data 
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NEAR FIELD DATA (Rec Ele #6) with k = 20 ; range = 14658 meters; array length = 1.885 meters 



Figure 117 - Near Field Data - computed and simulated 
Using 1.466 meter range, 1.571 meter long, 11 element near field array data 


FAR FIELD DATA (on axis) with k = 20; range = 1130.9734 meters; array length = 1.885 meters 



Figure 118 - Far Field Data - computed and simulated 
Using 1.466 meter range, 1.571 meter long, 11 element near field array data 
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Figure 119 - Ma-new and Ma (original) Data Sets 


Using 1.466 meter range, 1.571 meter long, 11 element near field array data 
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NEAR FIELD DATA (Rec Ele #8) with k = 20 ; range = 14658 meters; array length = 1.885 meters 



Figure 120 - Near Field Data - computed and simulated 
Using 1.466 meter range, 2.199 meter long, 15 element near field array data 


FAR FIELD DATA (on axis) with k = 20; range = 1130.9734 meters; array length = 1.885 meters 



Figure 121 - Far Field Data - computed and simulated 
Using 1.466 meter range, 2.199 meter long, 15 element near field array data 


116 
























Figure 122 - Ma-newand Ma (original) Data Sets 


Using 1.466 meter range, 2.199 meter long, 15 element near field array data 
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NEAR FIELD DATA (Rec Ele #7) with k = 20 : range = 1.885 meters; array length = 1.885 meters 



Figure 123 - Near Field Data - computed and simulated 
Using 1.885 meter range, 1.885 meter long, 13 element near field array data 


FAR FIELD DATA (on axis) with k = 20; range = 1130.9734 meters; array length = 1.885 meters 



Figure 124 - Far Field Data - computed and simulated 
Using 1.885 meter range, 1.885 meter long, 13 element near field array data 
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Figure 125 - Ma-newand Ma (original) Data Sets 
Using 1.885 meter range, 1.885 meter long, 13 element near field array data 
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APPENDIX A - EXAMPLE AB MATRIX 


The following pages show sections of a full augmented AB_Matrix as identified 
in this effort. This example is provided to assist in visualizing the matrix components, 
their relationship to one another and the overall circulant nature of the matrix. The 
example provided details an enclosing sphere of very small order. Due to this small size, 
this matrix would not be suited for actual use but makes it ideal for instructional 
purposes. 

For this example, an enclosing sphere was constructed with the following 
dimensions: 

numv = 5 this represents the number of patches in the “north-south” direction 

sph_sym =10 this represents the number of “slices” that comprise the sphere 

numh = 3 this represents the number of elements in the near field linear array 

The enclosing sphere so obtained would look like the following: 



Figure A1 - CHIEF Enclosing Sphere 
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As discussed in the write-up, the eonstituents of the AB Matrix are of the 
defined sizes listed below: 

Matrix 

A 
B 

delG 
G 

PE 
VE 

where: 

PE = Pressure equivalence matrix 
VE = Veloeity equivalenee matrix 

The AB_Matrix is formed from the eoneatenation of the above components. 
This matrix takes the form of: 


Ai 


A 2 

B 2 

Aso 

B 50 

delGi 

Gi 

delG 2 

G 2 

delGso Gso 

PEi 

0 

PE 2 

0 

PEso 

0 

0 

VEi 

0 

VE 2 

0 

VE 50 


The size of one “slice” of the AB_matrix is represented by the above 
arrangement. Eor the given example parameters, this matrix is of size 16 x 100. Since 
the full matrix is cireulant due to the symmetry of the ehosen enelosing surfaee (a 
sphere), the full AB_matrix is of size 160 x 100. (A matrix is considered to be eireulant 


Size (rows x columns) 

Example Matrix 

numv X (sph numv * sph sym) 

(4 it. 


5 x50 

5x50 

numh X (sph numv * sph sym) 


3 X 50 

3 X 50 

( sph numv-l) x (sph numv * sph sym) 

a a 

4x50 

4x50 
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if each of its rows, or columns, is a circular shift of the previous row, or column. If the 
matrix elements are considered to be matrices themselves and they obey the same 
circulant behavior as identified above, the matrix is more correctly identified as being 
block circulant.) 

The pages that follow show some sections of the full example AB_matrix. It was 
not practical nor of any real value to include the entire 16,000 element matrix. Rather, 
some representative subset was provided. The matrix has been numbered for reference 
purposes with column numbers appearing across the top or bottom and row numbers 
appearing on the left hand side. 

Since the array is circulant, the elements defined in the first slice (Rows 1 thru 16) 
and (Columns 1 thru 100) will repeat in a predictable manner throughout the rest of the 
matrix. Rows 1 thru 16 in Column 1 have been highlighted in red to show this repetition 
in subsequent slices. The data that appears in (Rows 1 thru 16), column 1 also appears in 
(Rows 17 thru 32), columnll and will appear in (Rows 33 thru 48), Column 21 and in 
other slices up until its final appearance in (Rows 145 thru 160), Column 91. 

Examination of the other matrix elements shows that subsequent slices are created 
by shifting the present slice sph sym (in our case this is 10) positions to the right and 
allowing those columns shifted past the last column (in our case this is 100) to be 
wrapped back to the beginning of the slice. This occurs a total of (sph_sym-l) times 
yielding the final matrix. 

Another significant point in the construction of this matrix is the manner in which 
the adjacent patch pressure equivalence (PE) and velocity equivalence (VE) is 
represented. Along any single row of PE or VE, there will appear a single 1 value 
followed immediately by a single -1 value; all other elements in these rows are zero. 
This forms the conditions of (pi - pi+i) or (vi - Vi+i) that is required for adjacent patch 
locations on the enclosing sphere. An examination of the matrix will actually show 1,0,- 
1 (as in row 9, column s 1 thru 3) that is used to implement this condition. The zero in 
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this case is actually associated with the veloeity term and thus does not interrupt the 
adjaceney requirement between sueeessive pressure terms. This is an artifact of the 
interleaving of the A and B matrices. 

The AB_matrix represents the known terms in the overdetermined system that it 
ereates. That is, in the example provided, there are 100 total unknowns. These 
eorrespond to 50 unknown pressures and 50 unknown veloeities. A single to-be-solved- 
for pressure that will multiply all elements of eolurnn 1 will be generated via the least 
squares solution teehnique. A single to-be-solved- for veloeity that will multiply all 
elements of column 2 will also be generated. Thus pressure values eorresponding to all 
odd numbered eolumns and veloeity values eorresponding to all even numbered eolumns 
are ultimately generated. 

Before leaving this appendix, it is signifieant to provide some additional 
eomments on the symmetry requirement of the aforementioned eomponents. In the 
MATLAB implementation, a single sliee of the augmented AB_matrix is formed. Its 
symmetry eonditions allow either the full AB_matrix to be determined (as was illustrated 
in this appendix) or it may be used direetly with “rotational solve” tools available in 
CHIEF for MATLAB to obtain a solution. This latter usage is preferred since it is 
eomputationally more effieient and imposes signifieantly less memory demands than the 
full defined AB_matrix. 

To utilize these “rotational solve” tools, however, requires that the AB_matrix 
sliee be representative of a seetion of the full bloek eireulant AB_matrix. It is for this 
latter reason that sueh eare was given to the definition of the enelosing sphere, the 
plaeement of the near field element points and the manner in which the adjaceney 
eonditions were defined. The symmetrical enclosing sphere satisfies the basie required 
eireulant behavior sinee it is related to the field points in a very predictable manner. That 
is, a rotation of the enclosing sphere by one sliee (with sph sym being the total number 
of increments in a eomplete rotation) yields the same values of A, B , G and delG with 
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the exception that their placement in the matrix has shifted. In the example of this 
appendix, the column s of data shift 10 columns to the right for each incremental rotation. 
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Augmented 
AB Matrix 


1 

2 

3 

4 

5 

6 

A and deIG 

B and G 

A and deIG 

B and G 

A and deIG 

B and G 


1 

-2.995 + 0.173761 

3.5684 + 0.92134i 

-0.093182 + 0.151741 

0.54411 -0.19015i 

-0.013249-0.358891 

-0.032079 + 0.62409i 

2 

-0.10782 + 0.00369071 

0.52906 + 0.11537i 

-2.6184 + 0.14508i 

1.0352 -0.82771 i 

-0.11272 + 0.20776i 

0.7104-0.440971 

3 

0.017724-0.102391 

-0.047656 + 0.21304i 

-0.12407 +0.16402i 

0.69458-0.36421 i 

-2.6275 - 0.083952i 

0.67693 + 0.59722i 

4 

0.023804 + 0.097903i 

-0.01765-0.1466i 

-0.02217-0.28557i 

0.00057447 + 0.51634i 

-0.11272 + 0.20776i 

0.7104-0.440971 

5 

0.17511 +0.27727i 

-0.17103-0.30889i 

-0.10198 + 0.0442711 

0.11917-0.098908i 

-0.013249-0.35889i 

-0.032079 + 0.62409i 

6 

4.6858 +2.4081 i 

-10.891 - 14.991i 

1.2014-6.6544i 

-4.6272 + 28.753i 

2.9742 - 6.3597i 

-26.025 + 12.142i 

7 

-3.0872 + 1.95491 

17.786 + 3.1931i 

4.5592 + 0.99307i 

-29.777 - 6.7274i 

6.4483-3.513i 

-29.578 - 5.4269i 

8 

-2.5756-0.0910331 

3.9504 + 17.378i 

1.8408+ 1.0087i 

-23.936- 17.67i 

2.9742 - 6.3597i 

-26.025 + 12.142i 

9 

1 

0 

-1 

0 

0 

0 

10 

0 

0 

1 

0 

-1 

0 

11 

0 

0 

0 

0 

1 

0 

12 

0 

0 

0 

0 

0 

0 

13 

0 

1 

0 

-1 

0 

0 

14 

0 

0 

0 

1 

0 

-1 

15 

0 

0 

0 

0 

0 

1 

16 

0 

0 

0 

0 

0 

0 

17 

0.12274-0.024667i 

-0.91857-0.719031 

0.12452 + 0.129671 

-0.15791 -0.45878i 

-0.25971 -0.19699i 

0.40766 + 0.35186i 

18 

-0.0027186-0.048651 i 

-0.039756 +0.083106i 

-0.1168-0.163941 

0.45809 + 0.83399i 

0.05291 -0.014402i 

-0.1824-0.0310351 

19 

-0.028833 + 0.041861 i 

0.079527 - 0.06373i 

0.13329-0.062714i 

-0.34442 + 0.075509i 

-0.11712 + 0.22948i 

0.70567-0.51815i 

20 

0.14422-0.0026827i 

-0.18576-0.0217151 

-0.12969 + 0.00516171 

0.20041 - 0.0069979i 

0.05291 -0.014402i 

-0.1824-0.0310351 

21 

0.25895 +0.24221 i 

-0.26248 - 0.27254i 

-0.033291 + 0.026667i 

0.0081437 - 0.066227i 

-0.25971 -0.19699i 

0.40766 + 0.35186i 

22 

-3.4698 + 0.78644i 

8.1228 + 3.2625i 

-7.4585 + 0.94841 i 

10.618-6.5525i 

-18.52 + 17.577i 

24.984-26.491 i 

23 

0.38449 - 2.2907i 

-5.7626 + 2.9136i 

-3.2486-5.7139i 

9.1501 +9.1485i 

-26.567 + 4.8747i 

37.351 - 9.0999i 

24 

1.2807- 1.1678i 

-3.6595 - 2.6948i 

-2.6672-5.02121 

6.7898+ 10.278i 

-18.52 + 17.577i 

24.984-26.491 i 

25 

0 

0 

0 

0 

0 

0 

26 

0 

0 

0 

0 

0 

0 

27 

0 

0 

0 

0 

0 

0 

28 

0 

0 

0 

0 

0 

0 

29 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 

31 

0 

0 

0 

0 

0 

0 

32 

0 

0 

0 

0 

0 

0 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


7 

A and deIG 


8 

B and G 


9 

A and deIG 


10 

B and G 


11 

A and deIG 


12 

B and G 


-0.10198 + 0.044271 i 

0.11917-0.0989081 

0.17511 + 0.27727i 

-0.02217-0.285571 

0.00057447 + 0.51634i 

0.023804 + 0.097903i 

-0.12407 + 0.164021 

0.69458-0.364211 

0.017724-0.102391 

-2.6184 + 0.145081 

1.0352-0.827711 

-0.10782 + 0.00369071 

-0.093182 + 0.151741 

0.54411 -0.19015i 

-2.995 + 0.173761 

1.8408 + 1.0087i 

-23.936- 17.67i 

-2.5756-0.0910331 

4.5592 + 0.99307i 

-29.777 - 6.7274i 

-3.0872 + 1.9549i 

1.2014-6.65441 

-4.6272 + 28.753i 

4.6858 + 2.40811 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

1 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

1 

0 

-0.033291 + 0.026667i 

0.0081437 - 0.0662271 

0.25895 + 0.242211 

-0.12969 + 0.00516171 

0.20041 - 0.0069979i 

0.14422-0.00268271 

0.13329-0.0627141 

-0.34442 + 0.075509i 

-0.028833 + 0.041861 i 

-0.1168-0.163941 

0.45809 + 0.83399i 

-0.0027186-0.0486511 

0.12452 + 0.129671 

-0.15791 -0.45878i 

0.12274- 0.0246671 

-2.6672-5.02121 

6.7898 + 10.278i 

1.2807- 1.16781 

-3.2486-5.71391 

9.1501 + 9.1485i 

0.38449 - 2.2907i 

-7.4585 + 0.948411 

10.618-6.55251 

-3.4698 + 0.78644i 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


-0.17103-0.308891 

0.12274-0.0246671 

-0.91857-0.719031 

-0.01765-0.14661 

-0.0027186-0.0486511 

-0.039756 + 0.0831061 

-0.047656 + 0.213041 

-0.028833 + 0.041861 i 

0.079527 - 0.06373i 

0.52906 + 0.115371 

0.14422-0.00268271 

-0.18576-0.0217151 

3.5684 + 0.921341 

0.25895 + 0.242211 

-0.26248 - 0.27254i 

3.9504 + 17.378i 

-1.4722- 1.147i 

0.75928 + 3.8456i 

17.786 + 3.19311 

1.7204-0.407311 

-4.0252 - 2.8299i 

-10.891 - 14.991i 

1.3635 + 0.954431 

1.3403-6.03111 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

-0.26248 - 0.27254i 

-2.995 + 0.173761 

3.5684 + 0.921341 

-0.18576-0.0217151 

-0.10782 + 0.00369071 

0.52906 + 0.115371 

0.079527 - 0.06373i 

0.017724-0.102391 

-0.047656 + 0.213041 

-0.039756 + 0.0831061 

0.023804 + 0.0979031 

-0.01765-0.14661 

-0.91857-0.719031 

0.17511 +0.277271 

-0.17103-0.308891 

-3.6595 - 2.6948i 

4.6858 + 2.40811 

-10.891 - 14.991i 

-5.7626 + 2.91361 

-3.0872 + 1.95491 

17.786 + 3.19311 

8.1228 + 3.26251 

-2.5756-0.0910331 

3.9504+ 17.378i 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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19 

A and deIG 


20 

B and G 


21 

A and deIG 


22 

B and G 


23 

A and deIG 


24 

B and G 


1 

0.25895 + 0.242211 

-0.26248 - 0.27254i 

-0.016032-0.0141341 

0.044255 + 0.34347i 

-0.081084-0.159981 

0.094108 + 0.397591 

2 

0.14422-0.00268271 

-0.18576-0.0217151 

0.10769 + 0.179371 

-0.16384-0.447221 

-0.054417 + 0.0983391 

0.14141 -0.23786i 

3 

-0.028833 + 0.0418611 

0.079527 - 0.06373i 

-0.013898-0.248111 

-0.016261 +0.36104i 

-0.074295-0.156541 

0.13711 + 0.22453i 

4 

-0.0027186-0.0486511 

-0.039756 + 0.0831061 

-0.17666-0.267451 

0.18622 + 0.325061 

0.16616 + 0.0882261 

-0.18519-0.116331 

5 

0.12274-0.0246671 

-0.91857-0.719031 

0.41548 + 0.0498141 

-0.43229 - 0.070006i 

-0.12447-0.0801011 

0.10459 + 0.124611 

6 

1.3635 + 0.954431 

1.3403-6.03111 

0.77887 + 1.2599i 

10.297 + 4.61521 

-16.185-2.33461 

-22.008-3.81571 

7 

1.7204-0.407311 

-4.0252 - 2.8299i 

-2.9368 + 0.28087i 

-11.836 + 3.8171 

1.0893- 16.6911 

3.0386- 18.0041 

8 

-1.4722- 1.147i 

0.75928 + 3.8456i 

-3.1626-3.73351 

-9.4768 - 9.6625i 

-1.4769- 17.695i 

-2.173- 14.666i 

9 

0 

0 

0 

0 

0 

0 

10 

0 

0 

0 

0 

0 

0 

11 

0 

0 

0 

0 

0 

0 

12 

0 

0 

0 

0 

0 

0 

13 

0 

0 

0 

0 

0 

0 

14 

0 

0 

0 

0 

0 

0 

15 

0 

0 

0 

0 

0 

0 

16 

0 

0 

0 

0 

0 

0 

17 

0.17511 +0.27727i 

-0.17103-0.308891 

0.12274-0.0246671 

-0.91857-0.719031 

0.12452 + 0.12967i 

-0.15791 -0.45878i 

18 

0.023804 + 0.097903i 

-0.01765-0.14661 

-0.0027186-0.0486511 

-0.039756 + 0.083106i 

-0.1168-0.163941 

0.45809 + 0.83399i 

19 

0.017724-0.102391 

-0.047656 + 0.213041 

-0.028833 + 0.041861 i 

0.079527 - 0.06373i 

0.13329-0.0627141 

-0.34442 + 0.075509i 

20 

-0.10782 + 0.00369071 

0.52906 + 0.11537i 

0.14422-0.00268271 

-0.18576-0.0217151 

-0.12969 + 0.00516171 

0.20041 - 0.0069979i 

21 

-2.995 + 0.173761 

3.5684 + 0.92134i 

0.25895 + 0.242211 

-0.26248 - 0.27254i 

-0.033291 + 0.026667i 

0.0081437 - 0.0662271 

22 

-2.5756-0.0910331 

3.9504+ 17.378i 

-1.4722- 1.147i 

0.75928 + 3.8456i 

2.8948-3.31811 

3.8416- 13.271 

23 

-3.0872 + 1.9549i 

17.786 + 3.1931i 

1.7204-0.407311 

-4.0252 - 2.8299i 

5.1313 + 2.27991 

16.787 + 3.8307i 

24 

4.6858 + 2.40811 

-10.891 - 14.991i 

1.3635 + 0.954431 

1.3403-6.03111 

7.4033 + 3.3046i 

19.042 + 9.3051 i 

25 

0 

0 

0 

0 

0 

0 

26 

0 

0 

0 

0 

0 

0 

27 

0 

0 

0 

0 

0 

0 

28 

-1 

0 

0 

0 

0 

0 

29 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 

31 

0 

0 

0 

0 

0 

0 

32 

0 

-1 

0 

0 

0 

0 


128 



127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 


0 

0 

0.25895 +0.24221 i 
0.14422-0.0026827i 
-0.028833+ 0.041861 i 
-0.0027186-0.048651 i 
0.12274-0.024667i 
1.2807- 1.1678i 
0.38449 - 2.2907i 
-3.4698 + 0.78644i 
0 
0 
0 
0 
0 
0 
0 
0 

0.41548 + 0.049814i 
-0.17666-0.26745i 
-0.013898-0.24811i 
0.10769+ 0.17937i 
-0.016032-0.014134i 
-2.8241 + 4.7513i 
2.6929 + 6.5929i 
5.4551 - 6.5507i 
0 
0 
0 
0 
0 
0 
0 
0 

79 

A and deIG 


0 

-1 

-0.26248 - 0.27254i 
-0.18576-0.021715i 
0.079527 - 0.06373i 
-0.039756 + 0.0831061 
-0.91857-0.719031 
-3.6595 - 2.6948i 
-5.7626+ 2.9136i 
8.1228 + 3.26251 
0 
0 
0 
0 
0 
0 
0 
0 

-0.43229 - 0.070006i 
0.18622 + 0.32506i 
-0.016261 + 0.36104i 
-0.16384-0.447221 
0.044255 + 0.34347i 
6.5666- 10.689i 
-5.0602- 11.867i 
-8.0569 + 10.1371 
0 
0 
0 
0 
0 
0 
0 
0 

80 

B and G 


0 

0 

-2.995 + 0.173761 
-0.10782 + 0.00369071 
0.017724-0.102391 
0.023804 + 0.097903i 
0.17511 +0.27727i 
4.6858 + 2.40811 
-3.0872 + 1.95491 
-2.5756 -0.0910331 
1 
0 
0 
0 
0 
0 
0 
0 

0.12274-0.0246671 
-0.0027186-0.0486511 
-0.028833 + 0.041861 i 
0.14422-0.0026827i 
0.25895 + 0.242211 
-3.4698 + 0.78644i 
0.38449 - 2.2907i 
1.2807- 1.1678i 
0 
0 
0 
0 
0 
0 
0 
0 

81 

A and deIG 


0 

0 

3.5684 + 0.921341 
0.52906 + 0.115371 
-0.047656 + 0.213041 
-0.01765-0.1466i 
-0.17103-0.30889i 
-10.891 - 14.991i 
17.786+ 3.1931i 
3.9504+ 17.378i 
0 
0 
0 
0 
1 
0 
0 
0 

-0.91857-0.71903i 
-0.039756 +0.083106i 
0.079527 - 0.06373i 
-0.18576-0.021715i 
-0.26248 - 0.27254i 
8.1228 + 3.26251 
-5.7626 +2.9136i 
-3.6595 - 2.6948i 
0 
0 
0 
0 
0 
0 
0 
0 

82 

B and G 


0 

0 

-0.093182 + 0.15174i 
-2.6184 + 0.14508i 
-0.12407 + 0.16402i 
-0.02217-0.285571 
-0.10198 + 0.044271 i 
1.2014 -6.6544i 
4.5592 + 0.99307i 
1.8408 + 1.0087i 
-1 
1 
0 
0 
0 
0 
0 
0 

0.12452 + 0.12967i 
-0.1168-0.163941 
0.13329-0.062714i 
-0.12969+ 0.0051617i 
-0.033291 + 0.026667i 
-7.4585 + 0.94841 i 
-3.2486-5.7139i 
-2.6672-5.0212i 
0 
0 
0 
0 
0 
0 
0 
0 
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A and deIG 


0 

0 

0.54411 -0.19015i 
1.0352-0.827711 
0.69458-0.364211 
0.00057447 + 0.516341 
0.11917-0.0989081 
-4.6272 + 28.753i 
-29.777 - 6.7274i 
-23.936- 17.67i 
0 
0 
0 
0 
-1 
1 
0 
0 

-0.15791 -0.45878i 
0.45809 + 0.83399i 
-0.34442 + 0.075509i 
0.20041 - 0.0069979i 
0.0081437-0.0662271 
10.618-6.55251 
9.1501 +9.1485i 
6.7898 + 10.278i 
0 
0 
0 
0 
0 
0 
0 
0 
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B and G 


129 


127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 


0 

0 

0.12274-0.024667i 
-0.0027186-0.048651 i 
-0.028833 + 0.041861 i 
0.14422-0.0026827i 
0.25895 +0.24221 i 
-1.4722- 1.147i 
1.7204-0.40731 i 
1.3635 + 0.95443i 
0 
0 
0 
0 
0 
0 
0 
0 

-2.995 + 0.173761 
-0.10782 + 0.00369071 
0.017724-0.102391 
0.023804 + 0.097903i 
0.17511 +0.27727i 
4.6858 +2.4081 i 
-3.0872 + 1.95491 
-2.5756-0.0910331 
1 
0 
0 
0 
0 
0 
0 
0 
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A and deIG 


0 

0 

-0.91857-0.719031 
-0.039756 +0.083106i 
0.079527 - 0.06373i 
-0.18576-0.021715i 
-0.26248 - 0.27254i 
0.75928 + 3.8456i 
-4.0252 - 2.8299i 
1.3403-6.0311i 
0 
0 
0 
0 
0 
0 
0 
0 

3.5684+ 0.92134i 
0.52906 + 0.11537i 
-0.047656 + 0.21304i 
-0.01765-0.1466i 
-0.17103-0.30889i 
-10.891 - 14.991i 
17.786 + 3.1931i 
3.9504+ 17.378i 
0 
0 
0 
0 
1 
0 
0 
0 
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B and G 


0 

0 

0.12452 + 0.129671 
-0.1168-0.163941 
0.13329-0.0627141 
-0.12969 + 0.00516171 
-0.033291 + 0.0266671 
2.8948-3.3181i 
5.1313 + 2.27991 
7.4033 + 3.3046i 
0 
0 
0 
0 
0 
0 
0 
0 

-0.093182 + 0.151741 
-2.6184 + 0.14508i 
-0.12407 + 0.16402i 
-0.02217-0.28557i 
-0.10198 + 0.044271 i 
1.2014-6.6544i 
4.5592 + 0.99307i 
1.8408+ 1.0087i 
-1 
1 
0 
0 
0 
0 
0 
0 
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A and deIG 


0 

0 

-0.15791 -0.45878i 
0.45809 + 0.83399i 
-0.34442 + 0.075509i 
0.20041 - 0.0069979i 
0.0081437 - 0.066227i 
3.8416- 13.27i 
16.787+ 3.8307i 
19.042 + 9.3051i 
0 
0 
0 
0 
0 
0 
0 
0 

0.54411 -0.19015i 
1.0352-0.82771 i 
0.69458-0.36421 i 
0.00057447 + 0.516341 
0.11917-0.098908i 
-4.6272 + 28.753i 
-29.777 - 6.7274i 
-23.936- 17.67i 
0 
0 
0 
0 
-1 
1 
0 
0 

94 

B and G 


0 

0 

-0.25971 -0.19699i 
0.05291 - 0.014402i 
-0.11712 + 0.22948i 
0.05291 - 0.014402i 
-0.25971 -0.19699i 
7.2473 +21.658i 
-8.2609 + 23.1111 
7.2473 + 21.6581 
0 
0 
0 
0 
0 
0 
0 
0 

-0.013249-0.35889i 
-0.11272 + 0.20776i 
-2.6275 - 0.0839521 
-0.11272 + 0.20776i 
-0.013249-0.35889i 
2.9742 - 6.3597i 
6.4483-3.513i 
2.9742 - 6.3597i 
0 
-1 
1 
0 
0 
0 
0 
0 
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A and deIG 


0 

0 

0.40766 + 0.35186i 
-0.1824-0.0310351 
0.70567-0.518151 
-0.1824-0.0310351 
0.40766 + 0.35186i 
18.329 + 37.397i 
-9.7121 +43.426i 
18.329 + 37.3971 
0 
0 
0 
0 
0 
0 
0 
0 

-0.032079 + 0.62409i 
0.7104-0.440971 
0.67693 + 0.59722i 
0.7104-0.440971 
-0.032079 + 0.62409i 
-26.025 + 12.142i 
-29.578 - 5.4269i 
-26.025 + 12.142i 
0 
0 
0 
0 
0 
-1 
1 
0 
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B and G 


130 


127 

0 

0 

0 

0 

128 

0 

0 

0 

0 

129 

-0.033291 + 0.026667i 

0.0081437 - 0.066227i 

0.25895 + 0.242211 

-0.26248 - 0.27254i 

130 

-0.12969 + 0.0051617i 

0.20041 - 0.0069979i 

0.14422-0.0026827i 

-0.18576-0.021715i 

131 

0.13329-0.062714i 

-0.34442 + 0.075509i 

-0.028833 + 0.041861 i 

0.079527 - 0.06373i 

132 

-0.1168-0.16394i 

0.45809 + 0.83399i 

-0.0027186-0.0486511 

-0.039756 + 0.0831061 

133 

0.12452 + 0.12967i 

-0.15791 -0.45878i 

0.12274-0.0246671 

-0.91857-0.71903i 

134 

7.4033 + 3.3046i 

19.042 + 9.3051 i 

1.3635 + 0.954431 

1.3403-6.0311i 

135 

5.1313 + 2.27991 

16.787 + 3.8307i 

1.7204-0.40731 i 

-4.0252 - 2.8299i 

136 

2.8948-3.31811 

3.8416- 13.27i 

-1.4722- 1.147i 

0.75928 + 3.8456i 

137 

0 

0 

0 

0 

138 

0 

0 

0 

0 

139 

0 

0 

0 

0 

140 

0 

0 

0 

0 

141 

0 

0 

0 

0 

142 

0 

0 

0 

0 

143 

0 

0 

0 

0 

144 

0 

0 

0 

0 

145 

-0.10198 + 0.044271 i 

0.11917-0.098908i 

0.17511 +0.27727i 

-0.17103-0.30889i 

146 

-0.02217-0.285571 

0.00057447 + 0.516341 

0.023804 + 0.097903i 

-0.01765-0.1466i 

147 

-0.12407 + 0.16402i 

0.69458-0.36421 i 

0.017724-0.10239i 

-0.047656 + 0.213041 

148 

-2.6184 + 0.14508i 

1.0352-0.82771 i 

-0.10782 + 0.00369071 

0.52906 + 0.115371 

149 

-0.093182 + 0.15174i 

0.54411 -0.19015i 

-2.995 + 0.173761 

3.5684 + 0.921341 

150 

1.8408 + 1.0087i 

-23.936- 17.67i 

-2.5756-0.091033i 

3.9504 + 17.378i 

151 

4.5592 + 0.99307i 

-29.777 - 6.7274i 

-3.0872 + 1.9549i 

17.786 + 3.1931i 

152 

1.2014-6.65441 

-4.6272 + 28.753i 

4.6858 +2.4081 i 

-10.891 - 14.991i 

153 

0 

0 

0 

0 

154 

0 

0 

0 

0 

155 

-1 

0 

0 

0 

156 

1 

0 

-1 

0 

157 

0 

0 

0 

0 

158 

0 

0 

0 

0 

159 

0 

-1 

0 

0 

160 

0 

1 

0 

-1 


97 

98 

99 

100 


A and deIG 

B and G 

A and deIG 

B and G 
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APPENDIX B - LEAST SQUARES SOLUTION TECHNIQUE 


The least squares teehnique allows an overdetermined system (“m” equations in 
“n” unknowns, with m>n) to be evaluated in an optimized manner. This optimized 
manner results in values of x that minimize the quantity |Ax - b| for the originating 
equation of Ax = b. 

The veetor x that minimizes |Ax - b| is the solution to the normal equations: 

A^ Ax = A^ b 

T IT rvi 

This veetor x = (A*ArA*b is the least squares solution to Ax =b ^ k 

This appendix seeks to provide the reader with a little more understanding of the 
meehanism of this teehnique. A relatively simple example of a sinusoidal wave 
eorrupted with a higher frequeney random noise is generated. The least squares solution 
teehnique is then used to determine a best fit “pure” sine wave to the input data provided. 
The sine wave so obtained minimizes the error per the aforementioned solution 
definitions. 

The least squares solution technique can best be envisioned as a means to provide 
a “best” solution to a long list of equations. Each equation has the same number and type 
of unknown constants. For the example chosen, a sine wave can be written in terms of 
both sine and cosine components. That is, this signal can be written as the sum of these 
two entities via: 

D sin(cot) + E cos(cot) (1) 

where: 

CO = radian frequency 
t = time 
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D and E are some defined constants 


In addition, the signal may have a DC component, so Equation (1) may be augmented to: 

D sin(cot) + E cos(cot) + E (2) 

where: 

E = some constant DC level 

Equation (2) can uniquely specify a given radian frequency sinusoidally varying 
signal by defining the values for the constants D, E and E. Phase delays can be 
introduced by the proper weighting of the D and E values. That is, if E=0, then the signal 
so defined has pure sine characteristics. If D =0, then the signal so defined has pure 
cosine characteristics. Non-zero values for both D and E yield sinusoidally varying 
signals with some defined phase shift. 

Under the assumptions that the radian frequency, co, is a known value and that 
time, t, can be incrementally adjusted (ti, t 2 , t 3 , U, etc), a set of equations defining the 
signal may be created. 

D sin(coti) + E cos(coti) + E (3) 

D sin(cot 2 ) + E cos(cot 2 ) + E 
D sin(cot 3 ) + E cos(cot 3 ) + E 
D sin(cot 4 ) + E cos(cot 4 ) + E 


The equations defined above consist of 3 unknowns (D, E and E) that will 
eventually need to be solved for. To assist in this solution, the equations specified in (3) 
above are equated to a given a set of measured data values, m. This leads to: 
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(4) 


D sin(coti) + E cos(coti) + F = mi 
D sin(cot2) + E cos(cot2) + F = m2 
D sin(cot3) + E cos(cot3) + F = m3 
D sin(cot4) + E cos(cot4) + F = 114 


The expressions above form the overdetermined set of equations that will be 
solved via least squares solution teehnique. The equations in ( 4 ) are rewritten in the form 
of matriees yielding: 


sin(coti) 

cos(coti) 

1 *[DEF] ‘ 

= mi 

( 5 ) 

sin(cot2) 

cos(cot2) 

1 

= m2 


sin(cot3) 

cos(cot3) 

1 

= m3 


sin(cot4) 

cos(cot4) 

1 

= m 4 



where D multiplies all values in column 1 , E multiplies all values in column 2 and F 
multiplies all values in column 3 of the first matrix. 
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The matrices in Equation (5) are now in the form of the defining equation Ax =b with: 


A = 

sin(coti) 

cos(coti) 


sin(cot 2 ) 

cos(cot 2 ) 


sin(cot 3 ) 

cos(cot 3 ) 


sin(cot 4 ) 

cos(cot 4 ) 


b = mi 

m2 
m3 
rru 


and 

x= [DEE]’ 

The solution may now proceed via the least squares solution technique detailed at 
the beginning of this appendix. 

Eor the example chosen, a 1 Hz sinusoidally varying signal is sampled at every 
1/256 seconds. The data corresponding to the first 10 rows (of the 1024 total) of the A 
matrix above was computed to be: 
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A 


1 

0.02454 

0.99970 

1 

2 

0.04907 

0.99880 

1 

3 

0.07357 

0.99729 

1 

4 

0.09802 

0.99518 

1 

5 

0.12241 

0.99248 

1 

6 

0.14673 

0.98918 

1 

7 

0.17096 

0.98528 

1 

8 

0.19509 

0.98079 

1 

9 

0.21910 

0.97570 

1 

10 

0.24298 

0.97003 

1 


The data corresponding to the first 10 rows of the b matrix were determined to be: 

b 

1 0.68685 

2 0.48398 

3 0.64596 

4 0.64235 

5 0.81405 

6 0.80746 

7 0.74193 

8 0.63214 

9 0.93579 

10 0.84604 


In like fashion as matrix A, vector b is composed of 1024 total rows. The data 
presented in vector b consists of the sinusoidally varying signal corrupted by some higher 
frequency noise. This corruption by noise explains why the values in b do not vary in a 
purely sinusoidal manner. 
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The solution, x, for this set-up was determined to be: 


X 


D 1.50520 

E 0.00903 

F 0.50071 


Reeall that D represents the sine eomponent, E represents the cosine component 
and F represents the DC offset of the measured signal. 

The chart that follows shows the original signal (blue crosses) along with the least 
squares solution for a best fit sinusoidally varying signal (red line). 



Figure B1 - Feast Squares Solution for Sinewave with Noise and DC Offset 


It is significant here to address the requirement to have an over determined system 
to utilize the method of least squares. The sinusoidal signal discussed in this section was 
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represented by a system of 1024 equations in 3 unknowns. This led, quite obviously, to 
an over determined system (sinee m=1024 was mueh greater than n =3). 

The thesis method required the ineorporation of theoretieal and physieal 
parameters to aehieve an over determined status. The example afforded by the 
AB_Matrix eonstruetion of Appendix A ean best illustrate this. 

Appendix A deseribed an AB_Matrix with the following eonstituents: 

numv = 5 ; this represents the number of patehes in the “north-south” direetion 

sph sym = 10 ; this represents the number of “sliees” that eomprise the sphere 
numh = 3 ; this represents the number of elements in the near field linear array 

The eorresponding system required that a solution be determined for 50 individual 
patehes (numv x sph sym). Eaeh pateh required a pressure and veloeity solution; this, 
therefore, lead to the need to solve for a total of 100 unknowns (50 pressures and 50 
veloeities). 


1) The use of the first Helmholtz relation, 0 =-Ap+Bv, provides a total of 50 
equations that define the physieal system. The 50 equations result from the 
faet that the AB_Matrix is eireulant and is based upon the number of vertieal 
patehes (numv) and the symmetry order (sph sym) of the enelosing sphere. 
These 50 equations are, however, insuffieient for solution via least squares 
sinee 100 unknowns must be solved for. This system ineorporating just this 
Helmholtz relation is, therefore, under determined. 

2) The use of near field data was essential to the overall effort sinee physieal 
eharaeteristies of the system needed to be quantified. The ineorporation of 
these physieal eharaeteristies was embodied in the seeond Helmholtz relation, 
pressure at a field point = delG p + G v. The measured near field element 
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data with the above relationship formed another set of defining equations for 
the system. 

The equations formed from this relationship eorresponded to the number of 
elements in the near field array (numh =3) and the symmetry order 
(sph_sym=10). Thus, for the example, 30 equations were assoeiated with 
this Helmholtz relation. 

These 30 equations are added to the initial 50 equations from the previous step 
yielding 80 total equations defining the system. These 80 equations, however, 
still ean not solve for the 100 unknowns; the system is still under determined. 

In this example, a near field array of 3 elements was seleeted. It would have 
been possible to define a near field array with more elements that eould have 
met the over determined requirement of the least squares methodology (when 
used with the previous Helmholtz relations) . For example, a 10 element 
array would have yielded 100 equations. These 100 plus the previous 50 
would have yielded a suffieient number for a least squares solution. 

For this relatively simple example, a 10 element array would not have 
presented any real difficulties to implement. For other systems, the number of 
near field elements to meet the over determined requirement could approach 
staggering proportions. 

Rather than achieving the additional system definitions via more near field 
elements, a physics based requirement for pressures and velocities on the 
patches was examined. 

3) The patch adjacency conditions of pressure and velocity, pi = pi+i and Vi = 
Vi+i, provided the physical relationships to achieve both additional defining 
equations for the system and to assist in forcing a realizable solution. Each 
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adjacency condition provided a total of (numv-1) * sph sym equations. For 
the example, this yielded an additional 40 equations for both pressure and 
velocity. 


Thus, the three components above provided 50, 30, and 80 (40 and 40) defining 
equations, respeetively, for the system. Thus, an over determined system was aehieved 
(since 160 equations now existed to solve for 100 unknowns). 
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For completeness and further investigation , the MATLAB ending used to define and 
solve for the sinusoidal signal (eorrupted by noise) in this appendix is provided. 


%funetion linear_least_squares_fit_to_data_set_SINEWAVE 

%elear all 
elose all 

% Create Data Set 
max_points = 1024; 

time_faetor = 1/256; % sampling period 

faetor =3 % noise faetor - high numbers represent low noise - OPERATOR INPUT 
radian_phase = 0; % phase shift of eomputed data set - OPERATOR INPUT 
offset = 0.5; %DC offset of eomputed data set - OPERATOR INPUT 
waveform magnitude = 1.5; %magnitude of sinewave - OPERATOR INPUT 


for integer =I:max_pomts; 
x_value(integer) = integer; 
random_number=rand; 

magnitude(mteger) = waveform_magnitude.*sm((mteger/256)*2*pi*l+radian_phase) + 
((random_number-0.5)/faetor) + offset; %+ 2*sm((4*integer/256)*pi); 

end 

max_magnitude =sqrt(max(magnitude. ^2)) 
magnitude2 = magnitude; 

% Plot Original Computed Data Set ineluding any noise eontributions 
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figure(l) 

plot (x_value*time_factor, magnitude) 


% Plot Normalized Data Set 
figure(2) 

plot (x_value*time_faotor, magnitude!) 
omega = 2*pi 


time =[time_faetor:time_faetor:max_pomts*time_faetor] ; 
data = magnitude!; 

% Plot Normalized Data Set (to allow for direet eomparison to proeessed 

% data later) 

figure(3) 

plot( time,data,'b+') 
xlabel('time (t)') 
ylabel('data value') 

time_size = size(time,2) 

% Set-up Matriees and Solve for Least Square Coeffieients 
onesmatrix = ones([ 1 timesize]); 

A = [sin(omega*time) ; cos(omega*time); ones matrix]'; 


B = data'; 
Atranspose = A'; 
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Result_l = A_transpose * A; 
Result_2 = A transpose * B; 
Result_final = Result_l \ Result_2; 


Cvalue = Result_final(l) 

Dvalue = Result_final(2) 

Evalue = Result_final(3) 

textl = ['Least Squares Fit Gives:'] 

text2 = [ num2str(C_value) ' sin(wt) + ' num2str(D_value) ' eos(wt) + ' num2str(E_value)] 

position_x = (max(time) + min(time))/8 
position_y = (max(data)) % + min(data))/2 

% Plot data using least squares derived eoeffieients 
hold on 

plot(time,C_value*sin(omega*time) + D_value*oos(omega*time) + E_value, 'r-') 
text(position_x, position y, text2, 'FontSize', 14) 
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APPENDIX C - EXPLANATION OF MATLAB CODING 


This appendix supplements Chapter IV - Steps of the Methodology. Detailed 
explanations of these steps as implemented via the three pieees of MATLAB eoding are 
provided herein. 

drake_stepl_*,m 

This step aeeomplishes the following: 

1) It ereates a simulated horizontal transmitting line array of some user- 
defined length. The elements in the transmitting array are equally spaeed and are 
spaeed at 1/6 lambda maximum. Lambda represents the wavelength of the 
frequeney of operation of the array. This transmitting array definition would be 
bypassed in the event that an aetual deviee was used for aeoustie evaluation. In 
this latter ease, the maximum dimension of the aeoustie deviee (its length for a 
line array) would be entered as input. 

2) A near field vertieal array is also simulated in his seetion of the eode. 
This is the array that is used in the near field of the aeoustie deviee under test to 
aequire response data. In the eoding that is provided, the user may define an array 
of arbitrary length made up of some defined number of equally spaeed elements. 
In aetual praetiee, the data entries here would eorrespond to the aetual 
eonstruetion details of the near field array being used. 

3) An estimate of the near field-far field boundary for the defined aeoustie 
deviee is ealeulated. 

4) The definition of range for the plaeement of the near field vertieal array 
for the near field measurements is made. A far field range is also defined whieh 
plaees a single reeeive hydrophone well in the deviee under test’s far field. In the 
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existing coding the near field range is defined as some percentage of the distance 
between the enclosing sphere radius (described later) and the near field-far field 
boundary (near field range = enclosing sphere radius plus factorA*(nf-ff boundary 
- enclosing sphere radius); the far field range is represented as 100 times this 
latter factor. 

5) Test points corresponding to the (x,y,z) positions for each of the near field 
array vertical array elements rotated at some defined angular increment about the 
acoustic device under test are created. These test points form a cylinder about the 
acoustic device under test. Test points for the far field simulation are also created. 
In latter case, since the far field will be calculated from the response of a single 
hydrophone, these test points form a circle about the device under test, ft is 
important to note here that the enclosing sphere symmetry order and the number 
of angular increments in this section are equal. 

In addition to the test points defined here, another set of test points are 
defined to allow an evaluation of the overall technique when dimensional 
uncertainty exists. The user has the ability to enter in values of the error 
associated with the range to the near field array, error in its vertical placement and 
errors in its angular orientation. In the implementation provided, these test points 
(with known errors) are associated with simulated data values only. 

6) This coding also calculates the complex pressures seen at each of the near 
field vertical array’s elements at each defined test point and at each of the single 
far field hydrophone’s test points. The near field test point response data will be 
used to simulate actual measurement data; the far field data so computed will be 
used later to compare to the output from the methodology. The calculations in 
this step are created by the MATLAB routine “prsvel”. This is a custom function 
call that computes complex response given acoustic parameters and defined 
source and field points. The complex response at a particular field point is 
computed as the contribution from each defined source point. The complex 
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responses so obtained are converted to dB values and normalized for later display. 
This section would not be executed if actual measurement data existed. 

7) The CHIEF enclosing sphere is created within this piece of coding. This 
spherical surface encloses the acoustic device under test with a known 
geometrical configuration. This sphere will be diced up into a mosaic of small 
patches; a single normal pressure and velocity will be assigned to each of these 
patches at a later point. The radius of this enclosing sphere is deduced from the 
size of the acoustic device under test and proximity to the location where the near 
field array will be placed. In like fashion as the placement of the near field array, 
the coding was implemented to allow a percentage of the separation distance 
between the maximum radial dimension of the device under test and the location 
of the near field array to be defined (this is factorB in the coding). 

The number of vertical patches and the enclosing sphere’s symmetry order 
are also defined in this section. The number of vertical patches is based upon the 
assumption that lambda/6 spacing is a reasonable maximum criterion to minimize 
pressure and velocity changes between adjacent patches (this is examined in more 
detail below). The symmetry order is based upon the number of vertical patches 
with the additional constraint that the symmetry should be a highly composite 
number (a number that may be represented by powers of 2 times 3's and 5's). This 
latter requirement pertains to the EFT methodology of solving the matrix 
mathematics and will not be addressed in detail here. 

A plot of a single section of the enclosing sphere is made to visualize the 
enclosing surface. 

The lambda/6 criterion of patch spacing was selected based upon an 
empirical analysis. This analysis looked at the complex response at a point in the 
far field. This response was formed from the averaged sum of some defined 
number of common magnitude vectors that were equally spaced (in an arc) but all 
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had different orientations. The resulting response in the far field was examined 
for varying veetor spaeings (ie; varying fraetional lambda spaeings). This 
response was eompared to that which would be achieved if all vectors were 
perfectly aligned with the far field point (ie; all vectors pointed in the same 
direction). This analysis showed that lambda/6 spacing yielded an error of only 
about 0.4 dB in the far field response. This error was deemed as satisfactory for 
the intended goals of the technique and was also expected to help minimize the 
number of patches required for implementation. 

The figures below support the above explanation. Three spaeings 
(lambda/6, lambda/12 and lambda/4) are shown in these figures. The 
corresponding difference in far field magnitude from the first case where all 
vectors pointed in the same direction are summarized in these figures. In these 
examples, 50 vectors are shown in each figure. 
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Figure Cl - Vector Spacing Summary and Effect on Far Field Pressure 

8) Using CHIEF tools, the coding solves for the A and B matrices (see 
CHIEF - Introduction for details on these matrices). Unlike previous versions of 
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CHIEF, this version allows a cessation of computations after these matrices are 
solved for. 

9) The greens function integral function call is next implemented for all 
defined near field test points at a single angular orientation along with all 
enclosing sphere patches (pssph). This function call results in defined values of G 
and delG for later computations. 

10) The Equality matrix is next established for the pressure and velocity 
adjacency conditions. This matrix consists has dimensions of (sph numv-l x 
sph_sym*sph_numv). This arises because there are sph_numv-l equality 
conditions imposed; each condition equates the pressure or normal velocity at a 
patch to its immediately adjacent patch. This matrix takes the general form of: 


\— 1 

\— 1 

1 

0 

0 

0 ... 

0 

0 

0 

1 

\— 1 

1 

0 

0 ... 

0 

0 

0 

0 

0 

0 

0 ... 

1 

\— 1 
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The matrix above represents a section (sph_numv-l x sph numv) of the overall 
Equality matrix. The remaining matrix components are all zeros. An Equality 
matrix for both pressure and velocity adjacency conditions will be provided later. 

11) The final operation in this piece of coding establishes a single maximum 
value for the A matrix and a single maximum value from all computed G/rhoc 
and delG values. These values will be used in the next step to normalize the 
elements of the combined AB_Matrix and M matrices. 

12) The significant data from this step is saved to file. 
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drake_step2_*,m 


This step takes the data obtained from the first step as its input and provides the following 
funetions: 

1) The weighting faetors alphaO and betaO are defined by the user. The 
alphaO term is assoeiated with the theoretieal Helmholtz Integral relations while the 
betaO term is assoeiated with the near field element data. (In practiee, this latter 
data set would be as aequired by the near field array; in this effort the near field 
data set was simulated.) Sinee these values ean be selected by the user, an 
emphasis can be placed (when the pressure and velocity values of the enclosing 
sphere patches are determined) on either the theoretical expectations or on the 
empirical data obtained. These weighting factors are used along with the values 
from step 11 in the first piece of coding to create the terms alpha and beta. 

2) The augmented Aa matrix is created by concatenating the A* alpha, 
delG*beta, and Equality (pressure and velocity) matrices; the augmented Ba matrix 
is created by concatenating the B*alpha/rhoc; G*beta/rhoc and Equality (pressure 
and velocity) matrices. These two matrices are then combined in such a way to 
create the AB_Matrix of form [AlBl A2B2A3B3...] 

3) The augmented Ma matrix is next created by concatenating zeros 
(expected outputs for the Helmholtz Integral relations), the near field array element 
data, and zeros (corresponding to the expected outputs for the (pi - pi+i) and (vi - 
Vi+i) terms). 


4) The coding then uses the “rotational_symmetry_solve” function call to 
create the “pv solution” matrix. This function call is specific to the CHIEE for 
MATLAB implementation and allows a problem with rotational symmetry to be 
solved. Ax = B. A is a (M x N) matrix, shorthand for MxrN circulant matrix where 
r is the order of the symmetry; b is a (M by 1) matrix of rhs values. The matrix so 
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obtained, pv_solution, contains individually catered values of pressure and velocity 
for each patch on the enclosing sphere (as determined by the underlying least 
squares solution technique contained within the “rotational symmetry solve” 
function call.) 

5) It is not just sufficient for the solved for pressures and velocities to yield 
values (when greens function integrations are used) that agree with those acquired 
by the near field array elements. Without taking into account the adjacency 
conditions and the Helmholtz Integral relations, the solved for pressures and 
velocities could be physically unrealizable. To assess the pv solution matrix, the 
function call “rotational symmetry multiply” is used. This is another function that 
is specific to the MATLAB for CHIEF implementation. This function solves a 
problem with rotational symmetry, Ax=b. A is a (M x rN) matrix where r is the 
order of the symmetry and b is a (rN by 1) matrix of rhs values. The output from 
this operation is a new matrix designated as the Ma-new matrix. 

6) The balance of the coding in this section allows the operator to look at the 
constituents of both the Ma and Ma-new matrices and assess their agreement. 
Perhaps the most significant display in this set pertains to the full image outputs of 
the Ma and Ma-new matrices. These plots contain a lot of information that allow a 
direct comparison of the two sets of receive array element outputs and allow a 
comparison in level of the resulting Helmholtz Integral relations and adjacency 
conditions within the Ma-new matrix. 

7) The computations in this section can be repeated with new weighting 
factors (alphaO and betaO) applied to try and arrive at a best agreement between the 
two M matrices. These computations can be reiterated without the need to re¬ 
create the data set from the first step. 
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8) The primary output of this coding is the pv solution matrix which will be 
used to re-establish pressure and velocities on the enclosing sphere for ultimate 
extrapolation into the far field. 

drake_step3_*,m 

1) This coding takes the pv solution matrix from the previous code and splits 
it into the pressure and velocity components corresponding to each of the 
enclosing sphere patches. Computations using greens functions integrations are 
performed to arrive at the G and delG values using the enclosing sphere geometry 
and the mid-point element of the near field array field test points and the defined 
far field point. These values of G and delG are used along with the newly 
defiend pressures and velocities to compute pressure at a field point, 

p delG + V G. 

2) Overplots of the original near field data and that obtained from this 
technique are displayed. In addition, overplots of the original simulated far field 
data set and that obtained via this technique are also displayed. In practice, this 
last overplot would not be made, instead just the computed far field plot would be 
displayed (since the response would not be known apriori). 
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%step #1 - near field to far field extrapolation techniqne 

close all 

V.el=25;V.az=25; % defines reference view for 3-D plots 

% factor for placement of near field array: sph rad -l-factorA(far_field - sph rad) 
factorA=0.08 

% factor for placement of enclosing sphere: height/2 -l-factorB(range_near_field - height/2) 
factorB=0.1305 

%define acoustic media parameters and k/frequency values of operation for 

%the transmitting array under test 

ss=1500;rho=1000;rhoc=ss*rho; 

k=20 

freq=k*ss/(2*pi) 
lambda = ss/freq 

%define parameters for vertical near field array 

rec array spacing = [-3*lambda:lambda/2:3*lambda]; % vertical spacing in meters of 
% receive array elements 

max length rec array = max(rec_array_spacing) - min(rec_array_spacing) 
numrecarrayelements = size(rec_array_spacing,2) 

%define parameters for horizontal transmitting array (device under test) 
height=6*lambda % this is the length of line array in meters 

num array elements = (round(height/((l/6)*lambda)))-l-l; %no. of elements comprising array 

% calculate far field criteria for defined line array for info purposes 
farfieldcomputedcriteria = (height^2)/(lambda); 
maxcriteria = farfieldcomputedcriteria 

%define parameters for enclosing sphere based upon wavelength of 
%transmitting array and its length 

sph rad = (-(l-factorB)*height/2 -factorA*factorB*max_criteria)/... 
(factorB-(factorA*factorB)-1) 

%base number of latitude lines on lambda/6 criterion 
sph_numv=ceil(pi* 6 * sph rad/(lambda)) 

%this gives a highly composite number that speeds up FFT process later 
sph_sym=2*nearest_small_factor(sph_numv) 

%define location of vertical near field array 

range near field = -(-factorA*max_criteria - sph_rad*(l-factorA)) 

%define range where far field is guaranteed 

range far field =100* round((max_criteria)/lambda)*lambda %simulate far field 

% allow operator to assess number of patches and other parameters before 
% proceeding 
operator_choice=0; 

operator choice = input('Do you want to use these values? l=yes 0=no ') 
if operator choice == 1 % proceed with computations for a "yes" condition 
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%define radian angles for which near field data will be required (the 
%number of angles must match the enclosing sphere symmetry) 
th=[0:l:sph_sym-l]' *(360/(sph_sym))* pi/180; 

%define some degree of positional uncertainty to test methodology 

range error = 0; % error in meters in range to rec array 

rec array z error = 0; % error in meters for vertical placement of rec array 

th error = 0; % error in radians of rec array placement 

%define all near field point locations 

saved = zeros(size(th,l)*num_rec_array_elements,3); 

for counter=l :num_rec_array_elements 

test_pts = [range_near_field*sin(th),range_near_field*cos(th),... 

zeros(length(th), 1 )+rec_array_spacing(counter)]; 
test_pts_error = [(range_near_field+range_error)*sin(th+th_error)... 
,(range_near_field+range_error)*cos(th+th_error),... 
zeros(length(th),l)+rec_array_spacing(counter)+rec_array_z_error]; 
All_nf_test_pts( 1 +(counter-1 )*size(th) :counter*size(th), 1:3) =test_pts; 
All_nf_test_pts_w_error( 1 +(counter-1 )*size(th):counter*size(th), 1:3)... 
=test_pts_error; %these test points include positional errors 


end 

%defme all far field point locations 

test_pts_far=[range_far_field*sin(th),range_far_field*cos(th),zeros(length(th),l)]; 

% define all source points of the transmitting line array 

source_pts =gen_all_xyz_combinations(0,linspace(-height/2,height/2,num_array_elements),0); 


% 

%surface definition 

%create enclosing sphere around array 

pssph=patch_defme('uv_sphere',[-pi/sph_sym,pi/sph_sym,0,pi],l,l,l,sph_numv,... 
[sph_rad,sph_rad,sph_rad],[0,0,0], [0,0,0]); 

% p 1 sph contains patch data for entire enclosing sphere 
p 1 sph=patches_by_rotation(pssph,sph_sym,'z'); 

% display single slice of the enclosing sphere 
M=plot3_patches(pssph,V); 

% define test points for the entire enclosing sphere based upon patch centers 
[test_pts_encl_sph,norm_encl_sphere,ar_encl_sphere] = patch_interp(plsph,.5,.5); 


% - 

%simulate pressures at the two defined ranges (near field and far field) 

%define normal vector for all defined source points 
normals = ones(size(source_pts)); 
normals(:,l) = 0.*normals(:,l); 
normals(:,3) = O.*normals(:,2); 

P_prsvel_near_saved = []; 

% compute pressures and velocities for all defined near field test points 
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for couiiter=l :num_rec_array_elements 

near_field_poiiits = All_nf_test_pts_w_error( 1 +(couiiter-1 )*size(th):... 
counter* size(th), 1:3); 

[P_prsvel_near,dummy]=prsvel(source_pts,near_field_poiiits,k,rhoc,normals); 
P_prsvel_near_saved( 1 +(counter-1 )*num_array_elements:counter*... 
num_array_elements,l :size(th))= P_prsvel_near; 


end 

% compute pressures and velocities for all defined far field test points 
[P_prsvel_far,dummy]=prsvel(source_pts,test_pts_far,k,rlioc,normals); 


P_prsvel_near_summed_saved = []; 
for counter=l :num_rec_array_elements 

summedvalue = sum(P_prsvel_near_saved( 1 +(counter-1) *num_array_elements:.. 

counter*num_array_elements, 1: size(tli, 1))); 
P_prsvel_near_summed_saved(counter,:) = summedvalue; 
end 

P_prsvel_far_summed = sum(P_prsvel_far); 

AlldBPreportednearall = []; 

for counter=l :num_rec_array_elements 

entry = P_prsvel_near_summed_saved(counter,:); 
All_dBP_reported_near_all(counter,:) = dB(entry) - dB(max(entry)); 
max_rec_array_values(counter,l) = dB(max(entry)); 
end 

dBPreportedfar = db(P_prsvel_far_summed)-db(max(P_prsvel_far_summed)); 

maxfarfieldvalue = db(max(P_prsvel_far_summed)); 

% - 


%Plot dB pressures (simulated) for all near field test points 

for counter=l :num_rec_array_elements 
figure 

plot(th*l 80/pi,All_dBP_reported_near_all(counter,:)); 
axis([0 360 -80 0]); 

text2 = ['NEAR FIELD DATA (Rec Ele #' num2str(counter)') with k = ' num2str(k) 
'; range = ' num2str(range_near_field)' meters; array length = '... 
num2str(height)' meters']; 

text3 = [' far field criteria = ' num2str(max_criteria)' meters; frequency = '... 
num2str(freq)' Hz; lambda = ' num2str(lambda)' meters' ]; 

title(text2, 'FontSize', 12) 
text(0, -10, text3, 'FontSize', 8) 
xlabel('angle (degrees)'); 
ylabel('pressure in db'); 
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end 


% Plot dB pressures (simulated) for all far field test points 
figure 

plot(th* 180/pi,dBPreportedfar); 
axis([0 360 -80 0]); 

text2 = ['FAR FIELD DATA with k = ' num2str(k)'; range = ' num2str(range_far_field)... 
' meters; array length = ' num2str(height)' meters']; 

title(text2, 'FontSize', 12); 
text(0, -10, texts, 'FontSize', 8); 
xlabel('angle (degrees)'); 
ylabel('pressure in db'); 

% set up CHIEF parameters to eompute A and B matriees 
elear O 

0.symmetry='ROTATE'; 

0.symmetry_order= {sph_sym,'z'}; 

O. quadra ture_order=2; 

0.adaptive_quadrature=adaptive_quadrature_strategy(pssph); 

0.rho=rho; 

0.e=ss; 

0.frequeney=freq; 

0.matriees_only= 1; 

[p,v]=solve_chief(pssph,0); 

global A_ B_ 

N=size(A_, 1) ;M=size(A_,2); 

% set up to eompute greens funetion integrals for all defined near field 

% test points 

All_g_nf=[]; 

All_dg_nf=[]; 

for eounter=l :num_ree_array_elements 

test_pts = All_nf_test_pts(H-(eounter-l)*size(th):eounter*size(th),l :3); 
[g_nf,dg_nf|=greens_fimetion_integral(pssph,0,test_pts(l,:)); 

All_g_nf(l:(sph_sym*sph_numv),counter) = g_nf; 
All_dg_nf(l:(sph_sym*sph_numv),eounter) = dgnf; 

end 

% Set-up Equality Matrix Strueture 

Ezeros = zeros(sph_numv-l,sph_numv); 

Eequabty = zeros(sph_numv-l,sph_numv); 

for eounter=l:(sph_numv-l) 

E_equality(eounter,eounter)=l; 

Eequabty (eounter,eounter-i-1 )=-1; 
end 

E temp = [E equabty E zeros]; 
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for counter = l:sph_sym-2 
E temp = [E temp E zeros]; 
end 

Efmal = Etemp; 

Efmalzeros = zeros(size(E_fmal)); 

% Determine maximum values in A matrix and in G and delG matrices 
factorA = max(max(abs(A_))) 

max_g_nf_templ= max(abs(All_g_nf(l :(sph_sym*sph_numv),l))); 
maxdgnftemp 1 =max(abs(All_dg_nf( 1: (sphsym* sphnumv), 1))); 

for counter=2: (num rec array elements) 

max_g_nf_temp2 = max(abs(All_g_nf(l :(sph_sym*sph_numv),counter))); 

max g nf temp 1 =max(max_g_nf_temp 1, max_g_nf_temp2); 

max_dg_nf_temp2= max(abs(All_dg_nf(l:(sph_sym*sph_numv),counter))); 
max dg nf templ =max(max_dg_nf_templ, max_dg_nf_temp2); 

end 

maxgnf = abs(max_g_nf_templ); 
maxdgnf = abs(max_dg_nf_templ); 

% use the following factor to normalize the G and delG data sets (note that 

% G also gets divided by rhoc 

factorgdgnf = max(max_dg_nf,max_g_nf/rhoc); 

%clean up variables that are no longer needed to save memory 

clear saved test_ptts source_pts near_field_points dummy summed value P_prsvel_near_saved... 
P_prsvel_far g_nf dg nf E temp entry E equality P_prsvel_near E zeros 

%save the data set to use in the next step 

save drake full data set 

end 
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%step #2 - near field to far field extrapolation techniqne 


close all 

% load data set from step #1 
load drake full data set 

% set to 0 for first run, set to any other number for subsequent runs of step2 to save decision data 
first_run=l; 

if first_run==0 

saved_step2_data=[]; 

else 

load step2_decision_data 
end 


%for group=l :size(th) % set group number for data to be displayed - group #1 
% is the first set of hydrohone data and its zeros range for group is 1 to size(th) 
group =1 

% Set alphaO and betaO factors (these will be iterated for final solution) 

%betaO places emphasis on the empirical near field data 
betaO = 1000 

% alphaO places emphasis on the theoretical Helmholtz Integral data 
alphaO = 1 

% Compute final sclaing factors that will be used 
alpha = alphaO/factorA 
beta = betaO/factorgdgnf 

%create the augmented A and B matrices 

new A single slice = [-A_*alpha;All_dg_nf'*beta; E final; E final zeros]; 
new B single slice = [B_*alpha/rhoc;All_g_nf'*beta/rhoc; E_final_zeros;E_final]; 

% interleave columns of new A and new B single slices into new AB single slice 
ABMatrix = [newAsingleslice newBsingleslice]; 
clear new A single slice new B single slice 

% now, put indices in order to keep circulant: 

% interleave from [A1A2A3...B1B2B3...] to [A1B1A2B2A3B3...] (block 
% circulant form), where each Am amd Bn is numv by numv in size, 
perm_ix=reshape([l :M;M+[1 :M]],2*M,1); 

AB_]VIatrix=AB_]VIatrix(: ,perm_ix); 


% Set up rhs of equation , Ma, the vector of actual measurements 

% Define two vectors of zeros corresponding to A and B (number of rows) and 
% the Equality Matrix (number of rows) 

ZeroO = zeros(sph_numv,l); 

Zerol = zeros(sph_numv-l,l); 

Mtemp = []; 

% build up M matrix using ZeroO, Zerol and the measured near field 
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% pressures at each receive element 


for counter =l:size(th); 

% "P_prsvel_near_summed" represents the complex pressure measured at a 
% particular field point. In this code, it was computed via "prsvel") 

M temp = [M temp; ZeroO;P_prsvel_near_summed_saved(:,counter); Zerol; Zerol]; 
end 

Mfinal = M_temp*beta; 
clear M temp 

group size = (3*sph_numv) -2 + size(rec_array_spacing,2); 

%solve for pressure and velocity 

pv solution = rotational_symmetry_solve(AB_Matrix, M final, sph sym); 

%compute Ma-new using solved for pressures and velocities 

M_from_pv_soln=rotational_symmetry_multiply(AB_Matrix,pv_solution,sph_sym); 

% separate constituents from the two M matrices 
M final resizel = reshape(M_final, group size, sph sym); 

M_final_resize2 = M_final_resizel(sph_numv+1 :sph_numv+size(rec_array_spacing,2),:); 

M_final_resize3 = reshape(M_final_resize2,(size(rec_array_spacing,2)*sph_sym), 1); 

M_final_zero0_resize2 = M_fmal_resizel(l:sph_numv,:); 

M_final_zero0_resize3 = reshape(M_final_zeroO_resize2,(sph_numv*sph_sym), 1); 

Mfinalzero 1 _resize2 = 

Mfmalresize 1 (sph_numv+size(rec_array_spacing,2)+l :3 *sph_numv+size(rec_array_spacing,2)-2,:); 
M_final_zerol_resize3 = reshape(M_final_zerol_resize2,((2*sph_numv-2)*sph_sym), 1); 

M_from_pv_soln_resizel = reshape(M_from_pv_soln, group size, sph sym); 

M_from_pv_soln_resize2 = 

M_from_pv_soln_resize 1 (sph_numv+1: sph_numv+size(rec_array_spacing,2),:); 

M_from_pv_soln_resize3 = reshape(M_from_pv_soln_resize2,(size(rec_array_spacing,2)*sph_sym), 1); 

M_from_pv_soln_zero0_resize2 = M_from_pv_soln_resizel(l :sph_numv,:); 

M_from_pv_soln_zero0_resize3 = reshape(M_from_pv_soln_zeroO_resize2,(sph_numv*sph_sym),l); 

M_from_pv_soln_zerol_resize2 = 

M_from_pv_soln_resizel(sph_numv+size(rec_array_spacing,2)+l:3*sph_numv+size(rec_array_spacing,2) 

-2,:); 

M_from_pv_soln_zerol_resize3 = reshape(M_from_pv_soln_zerol_resize2,((2*sph_numv- 

2)*sph_sym),l); 

maxhydro = 10*(round((max(dB(abs(M_fmal_resize3))))/10)) 

ymax = maxhydro; 
ymin = max hydro - 150; 

figure 

plot(l:size(rec_array_spacing,2)*sph_sym,dB(abs(M_final_resize3)),'k.-', 
l:size(rec_array_spacing,2)*sph_sym,dB(abs(M_from_pv_soln_resize3)),'r.-') 

textlO = ['ALL HYDROPHONE DATA (dB); alphaO = ' num2str(alpha0) ' ; alpha = ' num2str(alpha) ' ; 
betaO = ' num2str(beta0)'; beta = ' num2str(beta)]; 
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title(textlO) 

axis([l size(rec_array_spacing,2)*sph_sym ymin ymax]) 
xlabel('Near Field Array Field Points') 
ylabel('Level (dB)') 

figure 

plot(th* 180/pi,min(dB(abs(M_from_pv_soln_resize2))),'r.-') 

textlO = ['min(Ma-new FIYDROPFIONE DATA) (dB); alphaO = ' num2str(alplia0) ' ; alpha = 

num2str(alpha)'; betaO = ' num2str(beta0)'; beta = ' num2str(beta)]; 

title(textlO) 

axis([0 400 ymin ymax]) 
xlabel('Angle (degrees)') 
ylabel('Level (dB)') 

M hydro difference = dB((]Vt_from_pv_soln_resize3 - M_final_resize3)/max(max(]Vt_final_resize3))); 
M_hydro_difference2 = dB(]Vt_from_pv_soln_resize3) - dB(]Vt_final_resize3); 

figure 

plot(l:size(rec_array_spacing,2)*sph_sym,M_hydro_difference,'g.-'); 

textlO = ['FIYDROPFIONE DATA - dB[(computed-aetual)/max(actual)] ; alphaO = ' num2str(alpha0) ' 

betaO =' num2str(beta0)]; 

title(textlO) 

xlabel('Near Field Array Field Points') 
ylabel('Level (dB)') 


figure 

plot(l:size(rec_array_spacing,2)*sph_sym,M_hydro_difference2,'g.-'); 

textlO = ['FIYDROPFIONE DATA - dB[(computed-actual)] ; alphaO = ' num2str(alpha0) ' ; betaO = 

num2str(beta0)]; 

title(textlO) 

xlabel('Near Field Array Field Points') 
ylabel('Level (dB)') 


figure 

plot(l:sph_numv*sph_sym,dB(abs(M_from_pv_soln_zero0_resize3)),'b.-') 

textlO = ['ALL Ma-new ZeroO Data (Helmholtz Integral) (dB); alphaO = ' num2str(alpha0) ' ; betaO = 

num2str(beta0)]; 

title(textlO) 

axis([l sph_numv*sph_sym ymin ymax]) 
ylabel('Level (dB)') 

xlabel('Helmholtz Integral Conditions (sphnumv x sphsym)') 
figure 

plot( 1: (2 * sph_numv-2) * sph_sym,dB(abs(M_from_pv_soln_zero 1 _resize3)), 'k. -') 

textlO = ['ALL Ma-new Zerol DATA (Pressure/Velocity Equivalence) (dB); alphaO = ' num2str(alpha0)' 

betaO =' num2str(beta0)]; 

title(textlO) 

axis([l (2*sph_numv-2)*sph_sym ymin ymax]) 
ylabel('Level (dB)') 

xlabel('Pressure/Velocity Adjacencies ((2*sphnumv-2) x sphsym)') 
figure 

plot(th* 180/pi,max(dB(abs(Mfrom_pv_soln_zeroO_resize2))),'b.-') 
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textlO = ['max(Ma-new ZeroO Data)- (Helmholtz Integral) (dB); alphaO = ' num2str(alpha0) ' ; betaO = ' 

num2str(beta0)]; 

title(textlO) 

axis([0 400 ymin ymax]) 
xlabel('Angle (degrees)') 
ylabel('Level (dB)') 

figure 

plot(th* 180/pi,max(dB(abs(M_from_pv_soln_zero l_resize2))),'k.-') 

textlO = ['max(Ma-new Zerol DATA) - (Pressure/Velocity Equivalence) (dB); alphaO = ' num2sb(alpha0)' 

; betaO = ' num2str(beta0)]; 

title(textlO) 

axis([0 400 ymin ymax]) 
xlabel('Angle (degrees)') 
ylabel('Level (dB)') 

M hydro zeroO diff = dB(min(M_from_pv_soln_resize2)) - dB(max(M_from_pv_soln_zeroO_resize2)); 

M hydro zerol diff = dB(min(M_from_pv_soln_resize2)) - dB(max(M_from_pv_soln_zerol_resize2)); 

figure 

plot(th* 180/pi,M_hydro_zeroO_diff,'g.-'); 

textlO = ['Ma-new [ min(HYDROPHONE) - max(ZeroO) ] (dB) ; alphaO = ' num2str(alpha0) ' ; betaO = ' 

num2sb(beta0)]; 

title(textlO) 

xlabel('Angle (degrees)') 
ylabel('Level (dB)') 

figure 

plot(th* 180/pi,M_hydro_zero l_diff,'g.-'); 

textlO = ['Ma-new [ min(HYDROPHONE) - max(Zerol) ] (dB) ; alphaO = ' num2str(alpha0) ' ; betaO = ' 

num2sb(beta0)]; 

title(textlO) 

xlabel('Angle (degrees)') 
ylabel('Level (dB)') 


figure 

plot(th* 180/pi,min(dB(abs(M_from_pv_soln_resize2))),'r.-',... 
th* 180/pi,max(dB(abs(M_from_pv_soln_zero0_resize2))),'b. 
th*180/pi,max(dB(abs(M_from_pv_soln_zerol_resize2))),'k.-') 
textlO = ['Ma-new - min(HYDROPHONE), max(ZeroO), max(Zerol); alphaO = ' num2sfr(alpha0) ' ; alpha 
= ' num2str(alpha)'; betaO = ' num2str(beta0)'; beta = ' num2sb(beta)]; 
textll = [' red: min(hydrophone); blue: max(ZeroO); black: max(Zerol)']; 

title(textlO) 

axis([0 400 ymin ymax]) 
xlabel('Angle (degrees)') 
ylabel('Level (dB)') 
text(0, ymin-i-10, textll, 'FontSize', 8) 

dim =[ymin,ymax] 
figure 

subplot(l,2,l);imagesc(max(-200,db(M_from_pv_soln_resizel)),cbm) 

colorbar 

title('Ma-new - from pv solution') 
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subplot(l,2,2);imagesc(max(-200,db(M_final_resizel)),clim) 

colorbar 

title('Ma - original') 

%clean up variables that are no longer needed 

clear A_ B_ All dg nf All g nf E final E fmal zeros M final M_from_pv_soln 


% compute parameters to be saved and examined later 

savedl = min(min(dB(abs(M_from_pv_soln_resize2)))) %minimum of the new hydrophone data minimum 
set 

saved2 = max(min(dB(abs(M_from_pv_soln_resize2)))) % maximum of the new hydrophone data 
minimum set 

saved3 = mean(min(dB(abs(M_from_pv_soln_resize2)))) %average of the new hydrophone data minimum 
set 

saved4 = min(max(dB(abs(M_from_pv_soln_zeroO_resize2)))) %minimum of the new zeroO data 

maximized set 

savedS = max(max(dB(abs(M_from_pv_soln_zeroO_resize2)))) %maximum of the new zeroO data 
maximized set 

savedb = mean(max(dB(abs(M_from_pv_soln_zeroO_resize2)))) %average of the new zeroO data 

maximized set 

saved? = min(max(dB(abs(M_from_pv_soln_zerol_resize2)))) %minimum of the new zerol data 

maximized set 

savedS = max(max(dB(abs(M_from_pv_soln_zerol_resize2)))) %maximum of the new zerol data 
maximized set 

saved9 = mean(max(dB(abs(M_from_pv_soln_zerol_resize2)))) %average of the new zerol data 

maximized set 

savedlO = min(M_hydro_difference2) %minimum of the (new-original) hydrophone data set 
savedl 1 = max(M_hydro_difference2) %maximum of the (new-original) hydrophone data set 
savedl2 = mean(M_hydro_difference2) %average of the (new-original) hydrophone data set 

savedl3 = min(M_hydro_zeroO_diff); %minimum of the (hydrophone(new) - zeroO(new) data 
savedl4 = max(M_hydro_zeroO_diff); %maximum of the (hydrophone(new) - zeroO(new) data 

saved 15 = min(M_hydro_zerol_diff); %minimum of the (hydrophone(new) - zerol (new) data 
savedlb = max(M_hydro_zerol_diff); %maximum of the (hydrophone(new) - zerol(new) data 


% order of data - 

% alphaO (1) 

%betaO (2) 

% min hydro diff (3) 

% max hydro diff (4) 

% mean hydro diff (5) 

% min hydro minus zeroO (6) 
% max hydro minus zeroO (7) 
% min hydro minus zero 1 (8) 
% max hydro minus zero 1 (9) 
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% min new hydro set (10) 

% max_new_hydro_set (11) 

% mean_new_hydro_set (12) 

% min new zeroO set (13) 

% max_new_zeroO_set (14) 

% mean_new_zero0_set(15) 

% min new zerol set (16) 

% max_new_zerol_set (17) 

% mean_new_zerol_set (18) 

step2_data = [alphaO betaO saved 10 saved 11 saved 12 saved 13 saved 14 ... 
saved 15 saved 16 saved 1 saved2 saved3 saved4 saved5 saved6 ... 
saved? savedS saved9]; 

if first run ==0 

saved_step2_data = [step2_data] 
else 

saved_step2_data = [saved_step2_data; step2_data]; 
end 

header_info={'alphaO (l)VbetaO (2) Vmin hydro diff (3)', 'max hydro diff (4)Vmean_hydro_diff (5)',... 
'min hyd - zeroO (6)',' max hyd - zeroO (7)Vmin_hyd_-_zerol (8)' 'max hyd - zerol (9)'} 

header=(l:9); 

[header; saved_step2_data(:, 1:9)] 


save step2_decision_data saved_step2_data; 
wklwrite('step2_data_wkl file', saved_step2_data) 
save drake full data set2 
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%step #3 - near field to far field extrapolation techniqne 


close all 

%load data set from step #2 
load drake full data set2 


% pull pressure and velocity values from previous step for use in 
% green's function integrations 

pv_solnl= resliape(pv_solution,2,size(pv_solution,l)/2); 

penclsphere = pv_solnl(l,:).'; 
venclsphere = pv_solnl(2,:).'/rlioc; 

clear pv solnl pv solution 


%select near field array element in the middle 
counter = ceil(num_rec_array_elements/2) 

%evaluate green's function from all patches to all near field points and 
%calculate pressures for all test points for the middle recieve array element 

test_pts = All_nf_test_pts(l+(counter-l)*size(th,l):counter*size(th,l),l :3); 

[g_nf,dg_nf|=greens_fiinction_integral(pssph,0,test_pts); 

fp_temp=calculate_field_pressures(dg_nf,g_nf,p_encl_sphere,v_encl_sphere); 

max_field_point = max(fp_temp); 

Field_pressures(counter,l:size(th,l)) = dB(fp_temp) - max_rec_array_values(counter,l); 
clear fp temp g_nf dg nf 

%Plot dB Pressures for near field (original and that obtained via 
%technique) 

figure 

plot(th*l80/pi,All_dBP_reported_near_all(counter,l:size(th,l)),'b.-', ... 
th* 180/pi,Field_pressures(counter, 1: size(th, 1)), 'r.:'); 

axis([0 360 -80 0]); 

text2 = ['NEAR FIELD DATA (Rec Ele #' num2str(counter)') with k = '... 

num2str(k)'; range = ' num2sfr(range_near_field)' meters; array length = '... 
num2str(height)' meters']; 

text3 = [' far field criteria = ' num2sfr(max_criteria)' meters; frequency = '... 

num2str(freq)' FIz; lambda = ' num2sfr(lambda)' meters' ]; 
text4 = [' alpha = ' num2sfr(alpha)'; beta = ' num2str(beta)'; alphaO = '... 
num2str(alpha0)'; betaO = ' num2str(beta0) ] 

text6 = [' number of equally spaced receive array elements: ' num2str(num_rec_array_elements)... 

'; its length is ' num2sfr(max_length_rec_array)' meters']; 
texts = [' red (computed via method); blue (measured/simulated)']; 
text? = [' enclosing sphere radius = ' num2sfr(sph_rad)' meters; sph numv = '... 
num2str(sph_numv)'; sph sym = ' num2sfr(sph_sym)]; 
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title(text2, 'FontSize', 12) 
text(0, -10, texts, 'FontSize', 8) 
text(0, -70, text4, 'FontSize', 8) 
text(0, -60, texts, 'FontSize', 8) 
text(0, -65, text6, 'FontSize', 8) 
text(0, -75, text7, 'FontSize', 8) 

xlabel('angle (degrees)'); 
ylabel('pressure in db'); 


%Plot dB Pressures differenees for near field (original and that obtained via 
%teehnique) 

figure 

plot(th*180/pi,All_dBP_reported_near_all(eounter,l:size(th,l))- Field_pressures(eounter,l:size(th,l)),'g.-'); 

text2 = ['NF DATA Differenee (Ree Ele #' num2str(eounter)') with k = ' 
num2str(k)'; range = ' num2str(range_near_field)' meters']; 

title(text2, 'FontSize', 12) 

xlabel('angle (degrees)'); 
ylabel('db'); 


%evaluate green's funetion and ealeulate pressure for all far field poinfs 

[g_ff,dg_ff|=greens_funefion_integral(pssph,0,tesf_pfs_far); 

P_ff_gfl=ealeulafe_field_pressures(dg_ff,g_ff,p_enel_sphere,v_enel_sphere); 

%normalize far field pressures eompufed via Greens Funetion and eonvert to dB 
%P_est_reported_ffl = db(P_ff_gfl)-db(max(P_ff_gfl)); 

Pestreportedffl = db(P_ff_gfl)-max_far_field_value; 

elear g_ff dgff 


%Plot dB Pressures for far field (original and that obtained via 

%teehnique) 

figure 

plot(th* 180/pi,dBP_reported_far,'b.-',th* 180/pi,P_est_reported_ffl, 'r.:'); 
axis([0 360 -80 0]); 

text2 = ['FAR FIELD DATA (on axis) with k = ' num2str(k)'; range = '... 

num2str(range_far_field)' meters; array length = ' num2str(height)' meters']; 
texts = [' far field eriteria = ' num2str(max_eriteria)' meters; frequeney = '... 

num2str(freq)' Flz; lambda = ' num2str(lambda)' meters' ]; 
text4 = [' alpha = ' num2str(alpha)'; beta = ' num2str(beta)'; alphaO = '... 

num2str(alpha0)'; betaO = ' num2str(beta0) ] 
text5 = [' red (eomputed via method); blue (measured/simulated)']; 

text6 = [' number of equally spaeed reeeive array elements: ' num2str(num_ree_array_elements)... 

'; its length is ' num2str(max_length_ree_array)' meters']; 
text7 = [' enelosing sphere radius = ' num2str(sph_rad)' meters; sph numv = '... 
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num2str(sph_numv)'; sph sym = ' num2str(sph_sym)]; 


title(text2, 'FontSize', 12) 
text(0, -10, texts, 'FontSize', 8) 
text(0, -70, text4, 'FontSize', 8) 
text(0, -60, texts, 'FontSize', 8) 
text(0, -65, text6, 'FontSize', 8) 
text(0, -75, text7, 'FontSize', 8) 

xlabel('angle (degrees)'); 
ylabel('pressure in db'); 


%Plot dB Pressures differenees for far field (original and that obtained via 

%teehnique) 

figure 

plot(th*180/pi,dBP_reported_far-P_est_reported_ffl, 'g.-'); 

text2 = ['FF DATA Differenee(on axis) with k = ' num2str(k)'; range = ' 
num2str(range_far_field)' meters']; 

title(text2, 'FontSize', 12) 

xlabel('angle (degrees)'); 
ylabel('db'); 


% Re-plot All Flydrophone Data that eorresponds to the 
% data set presented. 

figure 

plot(l:size(ree_array_spaeing,2)*sph_sym,abs(M_final_resize3),... 
l:size(ree_array_spaeing,2)*sph_sym,abs(M_from_pv_soln_resize3)) 

textlO = ['ALL HYDROPHONE DATA ; alphaO = ' num2str(alpha0)'; alpha 
num2str(alpha)'; betaO = ' num2str(beta0)'; beta = ' num2str(beta)]; 
title(textlO) 


elim =[ymin,ymax] 
figure 

subplot(l,2,l);imagese(max(-200,db(M_from_pv_soln_resizel)),ehm) 

eolorbar 

title('Ma-new - from pv solution') 

subplot(l,2,2);imagese(max(-200,db(M_final_resizel)),ehm) 

eolorbar 

title('Ma - original') 

range error % error in meters 
ree array z error % error in meters 
th error % error in radians 

save drake full data set final 
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